EVANS FUNCTIONS, JOST FUNCTIONS, 
AND FREDHOLM DETERMINANTS 



FRITZ GESZTESY, YURI LATUSHKIN, AND KONSTANTIN A. MAKAROV 

Abstract. The principal results of this paper consist of an intrinsic definition 
of the Evans function in terms of newly introduced generalized matrix- valued 
Jost solutions for general first-order matrix- valued differential equations on the 
real line, and a proof of the fact that the Evans function, a finite-dimensional 
determinant by construction, coincides with a modified Fredholm determinant 
associated with a Birman— Schwingcr-typo integral operator up to a nonvan- 
ishing factor. 



1. Introduction 

In this paper we study connections between the asymptotic behavior of solutions 
of the following first-order dxd matrix- valued systems of linear differential equations 
on the real line M, 

y'{x)^A{x)y{x), x€R, (1.1) 

y'ix) = iA{x) + R{x))y{x), a; S R, (1.2) 
which will often be referred to as the unperturbed and perturbed equation in the 
following. 

Three types of results are obtained: 

• First, we establish the existence of special matrix- valued solutions of the per- 
turbed equation (1.2), called the generalized matrix-valued Jost solutions, that are 
asymptotic to some reference solutions of the unperturbed equation (1.1). In spirit, 
these results resemble the celebrated Levinson Theorem [13, Theorem 1.3.1]; how- 
ever, we are using a different approach involving Bohl and Lyapunov exponents. 

• Second, we use the determinant of a finite-dimensional matrix composed of 
initial data of the generalized matrix- valued Jost solutions to calculate an infinite- 
dimensional modified Fredholm determinant of a Birman-Schwinger-type integral 
operator associated with equations (1.1) and (1.2). These results generalize a clas- 
sical relation identifying the Jost function and a Fredholm determinant in the case 
of one-dimensional half-line Schrodinger operators first derived by Jost and Pais 
[36]. 

• Third, using the generalized matrix-valued Jost solutions, we give a coordinate- 
free definition of the Evans function, and relate the Evans function to the above- 
mentioned finite-dimensional determinant. Simultaneously, this identifies the Evans 
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function with an infinite-dimensional modified Fredholm determinant associated 
with a Birman-Schwinger-type integral operator associated with (1.1) and (1.2) 
up to a nonvanishing factor. As a result, the Evans function is expressed via an 
infinite-dimensional modified Fredholm determinant, and, for the special case of the 
Schrodingcr equation, is proved to be equal to the classical Jost function familiar 
in scattering theory. 

The Evans function, E, is a Wronskian-type analytic function which is widely 
used to trace the spectrum of ordinary differential operators, £. Most frequently, 
the operators £ appear as linearizations of partial differential equations along spe- 
cial solutions such as travelling waves and steady states. For instance, C can be 
the one-dimensional Schrodinger operator, cP/dx"^ + V, obtained by linearizing the 
reaction-diffusion equation, Ut = Uxx + '''{u), about a steady state uq, so that 
V{x) = r'(wo(x)), X e K. 

Recently, the Evans function became one of the most important tools in stability 
analysis. It was originally introduced by J. W. Evans [15]-[18] to treat a particular 
model of nerve impulses, however, quite soon this object was generalized, and 
numerous connections of the Evans function to many fields of mathematical physics 
have been discovered (from Chern numbers in topology [1] to scattering data in 
quantum mechanics [39]). To give the reader just a small sample of the work done 
in this rapidly developing area, we cite, for instance, [1], [2], [4], [5], [11], [19], [20], 
[21], [22], [23]', [24], [32], [34], [35], [37], [38], [39]. [40], [42], [43], [47], [50], [51], [52], 
[54], [58], [59] , [60] , [61] , [62] , [67] , [68] , [69], and [70]. In addition, excellent reviews 
of this subject and further references can be found in [35] and [57]. 

Re-writing the eigenvalue problem for £ as a first-order system of differential 
equations, 

y'{x)=B{z,x)y{x), x e R, (1.3) 

where z G C is the spectral parameter, one is interested in conditions under which 
this system has C'-valued solutions y on M exponentially decaying at ±oo. The 
Evans function, E = E{z), is defined in such a way that E{z) = if and only if 
decaying C'-valucd solutions y exist (thus detecting the eigenvalues of C). Since 
in this paper we are not concerned with function theoretic properties of the Evans 
function such as its analyticity, etc., deferring this topic to a forthcoming pub- 
lication, wc simply fix a value zq of z and study equations (1.1) and (1.2) (and 
hence (1.3) with B{zo,x) = A{x) + R{x), x G M). We note that, according to all 
definitions available in the literature, the Evans function for (1.1) and (1.2) is not 
uniquely defined, and one of the central issues in the current paper is to motivate 
a "canonical" choice of the definition of E. 

To give an informal outline^ of our results, let us temporarily assume for sim- 
plicity of exposition, that A and R belong to L°°(R)'^^'^, and, in addition, that 
R € L^iR^'^. On the Hilbert space L'^{Ry consider operators, Ga and Ga+r, 
with domains i?^(M)'*, the usual Sobolev space, defined by 

{Gau){x) = -u'{x) + A{x)u{x), X e K, 

{Ga+ru){x) = -u'{x) + {A{x) + R{x))u{x), x G R. ^^''^^ 

Next, we represent the perturbation R as a product 

R{x) = Re{x)Rr{x), x G M. (1.5) 



See the glossary of notation at the end of this introduction. 
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Assuming the invertibihty of Ga and Ga+r, one can write 

Ga+r = Gj^ - G / RiRrG ji^^^, (1.6) 

and, in addition, one can show that the operator / + RrG^Ri is invertible. Mul- 
tiplying (1.6) from the left with R^. and solving for then yields 

G j^j^jf^ = G — G R(\I + RrG Ri] ^RrGj^. (l-^) 

Conversely, assuming the invertibihty of Ga and / + RrG^^Re, we conclude that 
Ga+r is invertible and (1.7) holds. Next, still assuming that Ga is invertible, one 
computes 

1/2^ ^—1/2 ^—1/2/^ r> D \y^^l/2 

A ^A+R(~'a —'^A a + niUr )Lr 

-1/2 -1/2 (-^-^^ 

= I + Gy^ RfRrGj^ , 

and assuming that R^G^^^ and RrG~^^^ are operators in Bi{L'^{S.Y) (for this it 
suffices that Rt,Rr € L'^{R)^'"'', cf. Lemma 2.9; also, see [30, 31, 66] for general 
information about the Schattcn von Neumann ideals Sp(-)), one obtains for the 
symmetrized perturbation determinant 

det2{G~^^^GA+RG^'/^) = dct2(/ + G-//^R,RrG^'^^) 

= det2{I + RrG^^Re). 

Here det2 (/ + T) denotes the modified Predholm determinant of a Hilbert-Schmidt 
operator T on L'^{R)'^ (and we used tlic fact that dct2(/ + TiTj) = det2(/ + T2T1) 
for Ti,T2 bounded operators such that T1T2 and T2T1 are Hilbert-Schmidt). In 
removing the restriction A,R£ Z/°°(]R)'^^'^ and permitting A to be locally inte- 
grable on R one needs to supply proper operator closures in formulas (1.6)-(1.9) 
(see, e.g., [27] and the extensive literature therein for a more general discussion of 
factorizations of perturbations). 

To make the connection with (1.1) and (1.2) we recall that the operator Ga 
is invertible on L^(M)'^ if and only if (1.1) has an exponential dichotomy Q on M 
(see, e.g., [7, Ch. 3]). By a classical result in the theory of dichotomic differential 
equations (see, e.g., [8, Proposition 4.1], [10, Ch. IV]), since (1.1) has an exponential 
dichotomy on M and the perturbation satisfies ||i?||c<ix<i e i^(R), one infers that 
the perturbed equation (1.2) has exponential dichotomies on both half-lines 1R+ 
and ]R_, and moreover (by yet another well-known Dichotomy Theorem, see [3], 
[48], [49], [56]), the operator Ga+r is then a Fredholm operator with zero Fredholm 
index. Thus, Ga+r is either invertible on L'^{M.y, or has a nontrivial null space 
(consisting exactly of the solutions of (1.2) that belong to H^(M.)'^). 

Due to the assumption ||i?||c<ixd G ^^(IR), the Birman-Schwinger-type operator 

K = RrG^^Ri (1.10) 

is a Hilbert-Schmidt operator on L^(M)'' (see Lemma 2.9). Thus, employing the 
celebrated Birman-Schwinger argument we conclude from (1.7) that the operator 
Ga+r is not invertible (that is, (1.2) has solutions in H^(M.)'^) if and only if —1 is 
an eigenvalue of K, or, in other words, if and only if the 2-modified perturbation 
determinant det2(/ -|- K) is equal to zero. 

Without going into further details at this point, wc mention that the Birman- 
Schwinger-type argument (cf. [55, Sect. X.III.3] for this terminology) in our present 
context provides a bijection between the geometric eigenspace (nuUspace) of Ga+r 
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and the geometric eigenspace of RrG^ Rg corresponding to the eigenvalue —1 (see, 
e.g., [25], [27], [55, Sect. XIII.3], [63, Ch. Ill], [64] and the recent discussion in [27] 
which particularly targets non- self- adjoint operators). 

Equation (1.9) finally illustrates the sense in which det2(/ + K) can be viewed 
as a properly symmetrized modified perturbation determinant (cf. [31, Sect. IV. 3]) 
associated with the pair (Ga, Ga+r)- We note at this point that det2(/ + K) will 
be one of the central objects in this paper. 

Our next step is to reduce the calculation of the infinite- dimensional determinant 
det2(/ -\- K) to & finite- dimensional determinant. This is possible because the 
integral kernel of the operator K happens to be of a special type: it belongs to 
the class of so-called semi-separable integral kernels (cf. [28], [29], [30, Ch. XIII]). 
For the integral operator with semi-separable integral kernel we construct a finite- 
dimensional determinant, D, such that the following representation holds: 

det2{I + K)=e®D, (1.11) 

where 8 is a certain constant. It is known from the work in [28] that D can be ex- 
pressed in terms of solutions of certain Volterra integral equations associated with 
(1.1) and (1.2). We make yet another step, and express D in terms of solutions 
of (1.2) that are asymptotic to some reference solutions of (1.1). We called these 
solutions of (1.2) the generalized matrix-valued Jost solutions, since they appear 
as a generalization of the classical Jost solutions of the Schrodinger equation [6, 
Ch. XVII], [46, Sect. 12.1]. In particular, we prove that the determinant D in 
(1.11) is equal to the classical Jost fimction in the case of the perturbed equa- 
tion (1.2) corresponding to the eigenvalue problem for the Schrodinger operator £. 
The reference solutions are chosen using the Lyapunov exponents while the rate 
of approximation of the reference solutions by the generalized matrix- valued Jost 
solutions is controlled by the Bohl exponents for (1.1) (see the terminology in [10, 
Ch. III]). The existence of the generalized matrix- valued Jost solutions is proved by 
a contraction mapping argument applied to a "mixed" system of integral equations 
of Volterra- and Fredholm-type (which appears to be a new element in the litera- 
ture). In addition, we identify the abstract asymptotic properties of the generalized 
matrix- valued Jost solutions leading to formula (1.11). Finally, wc complete the 
picture by proving that the finite-dimensional determinant D in (1.11) is equal to 
the Evans function when E is defined by means of the initial data of the generalized 
matrix-valued Jost solutions, thus making this choice of E canonical. 

We emphasize that our results on the existence of the generalized matrix- valued 
Jost solutions can be viewed (cf. Remarks 7.11 and 8.5) as a further refinement of 
the celebrated Levinson theorem (sec, for instance. [13, Theorem 1.3.1]). 

Some results of the current paper have been announced in [26] . A version of this 
theory for difference equations, based on the material developed in this paper, has 
been derived in [9]. 

The paper is organized as follows. In Section 2 we recall some known facts from 
[8, 10] regarding exponential dichotomies and the Bohl and Lyapunov exponents 
(which are our means to measure the asymptotic behavior of solutions of differ- 
ential equations), and also show that K is a, Hilbert-Schmidt operator on L^(M)'' 
(Lemma 2.9). In addition, we formulate Lemmas 2.13 and 2.14 (proved in Appendix 
B) to the effect that the exponential dichotomy on half-lines and, more generally, 
exponential splittings and the corresponding Bohl exponents for (1.1) persist un- 
der i^-perturbations. Note that the robustness of the Bohl exponents can also be 
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proved under assumptions different from ||_R||c<ixd £ L^(K) (e.g., assuming that the 
perturbation is continuous and decays at ±oc to zero in norm), and thus our subse- 
quent results can be developed in this different setting, but to keep this manuscript 
at a reasonable length we do not pursue this in the current paper. 

In Section 3 we introduce in Definition 3.1 the matrix- valued Jost solutions, Y± 
on R±, which are proved to exist (and to be unique) under the assumptions that 
(1.1) has an exponential dichotomy on R and ||i?||c<ix<i decays exponentially at ±oo 
at a rate controlled by the width of the two Bohl segments corresponding to the 
dichotomy projections (cf. Theorems 3.2 and 3.8). Although this result is obtained 
under strong exponential decay assumptions, it is shown to be optimal if one does 
not allow further exponential splitting of the dichotomy projections (cf. Example 
3.6). 

In Section 4 (cf. Theorem 4.1) we establish formula (1-11) with 

£> = detc.(> + (0) + r_(0)) (1.12) 

under the assumptions of Section 3. We note that these assumptions are satis- 
fied for compactly supported perturbations, and our method of proof of (1.11) is 
to establish this formula for a truncated perturbation first, and then pass to the 
limit as the support of the truncation expands. Formula (1.11) in the compactly 
supported perturbation case (which is our main intermediate calculation) is proved 
using results from [28] collected in Appendix A. We call D in (1.11) the Evans 
determinant. 

Section 5 deals with the class of perturbed equations (1.2), where each of the Bohl 
segments for the dichotomy projections degenerates to a single point (cf. Theorem 

5.3) . In particular, this assumption covers the case of (1.2) corresponding to the 
Schrodinger equation, and in this section we also show that our matrix-valued 
Jost solutions indeed generalize the classical Jost solutions, and, moreover, that 
the Evans determinant D in (1.11) coincides with the Jost function (see Theorem 

5.4) . We refer to [39, 40] for results relating the Evans function and a Fredholm 
determinant in the Schrodinger equation case. 

In Sections 6 and 7 we introduce and study the generalized matrix-valued Jost 
solutions, assuming that the exponential dichotomy for (1.1) admits further (and 
finer) exponential splitting (of order two in Section 6 and of an arbitrary order in 
Section 7). The generalized matrix- valued Jost solutions are defined abstractly as 
solutions of (1.2) that approximate decaying reference solutions of (1.1) at ±oo at 
an exponential rate which is better than the rate of decay of the reference solutions 
themselves (cf. Definitions 6.2 and 7.2). Unlike the matrix- valued Jost solutions 
from Section 3, the generalized matrix- valued Jost solutions arc not unique due to 
the possible appearance of "lower-order modes" , that is, solutions of the perturbed 
equation (1.2) having worse exponential rate of decay than the corresponding ref- 
erence sohitions. However, they arc unique up to terms that decay exponentially 
faster than the respective reference solutions (see Corollary 6.6 and Theorem 7.3). 
Nevertheless, we prove in Lemma 7.6 that the Evans determinant D defined by 
means of the generalized matrix- valued Jost solutions is determined uniquely. The 
existence of the generalized matrix- valued Jost solutions is proved using a "mixed" 
system of integral equations of Volterra- and Fredholm- type in Theorems 6.5 and 
7.3. The use of the finer exponential splitting allows us to prove the existence of 
the generalized matrix-valued Jost solutions and formula (1-11) under much weaker 
assumptions on the exponential decay of the perturbation than those in Section 3: 
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In Sections 6 and 7 the decay of the perturbation is controlled by the width of each 
of the Bohl segments corresponding to the fine exponential splitting (cf. Theorems 
6.5 and 7.8). This constitutes one of the main new effects observed in this paper: 
Passing from the case of two- to the case of many Bohl segments in the unperturbed 
equation (1.1) requires replacing the system of two Volterra integral equations by 
a cascade-type mixed system of Volterra- and Fredholm-type integral equations 
so that one can take care of the appearance of the above mentioned "lower-order 
modes". As a generalization in yet another direction, we express in Corollary 7.10 
the infinite-dimensional determinant det2(/ + K) in (1.11) as a ratio of two Evans 
determinants corresponding to the perturbed and the unperturbed equation (1.1) 
thus making the connection with the interpretation of det2(/ -|- K) as a (modified) 
perturbation determinant (cf. [31, Sect. IV. 3]). 

In Section 8 we treat the special but important case of the constant coefficient 
unperturbed equation (1.1), A{x) = A, x E M. We introduce the generalized 
matrix-valued Jost solutions as solutions of the "mixed" system of Volterra- and 
Predholm-type integral equations in Definition 8.2. Using the Jordan normal form of 
A, the existence of the solutions and formula (1.11) are proved under the following 
mild assumptions in Theorem 8.3: Either we suppose ||i?||ctixd G L^(IR) if the 
spectrum of A is semi-simple, or we assume a polynomial decay of ||i?||cdxd. These 
hypotheses arc shown to be optimal. 

Finally, in Section 9 we show that the Evans determinant D equals the Evans 
function for (1.1) and (1.2) defined via the generalized matrix-valued Jost solutions. 

A glossary of notation: For d,di,d2 G N, let C'^^^'^^ be the set of di x d2 
matrices with complex entries; by trcd(M) and detc<i(M) we denote the trace and 
determinant of a, d x d matrix M e C'*^''. || • He* denotes a vector norm in C^; 
II • llc-ix-i a matrix norm in C'*^''; || • H^- a norm in a Banach space X. For p > 1 
and J C R, LP{J), LP{JY, and LP{JY^'^ are the usual Lebesgue spaces on J with 
values in C, C'', and C''^'', associated with Lebesgue measure dx on J. Similarly, 
LP{J; w{x)dx), LP{J; w{x)dxY and L'^( J\ w(x)dxY^''' are the weighted spaces with 
a weight w > 0; the spaces of bounded continuous functions on J are denoted by 
Cb( J), Cb(J)'^ and Ch{JY'^'^, respectively. The identity matrix in C^^"* is denoted 
by Id and the identity operator on a Banach space X is denoted by / (or by Ix 
if its dependence on X needs to be stressed). For a projection Q we write ran((5) 
and ker((5) to denote the range and the null space (kernel) of Q. If T is a linear 
operator on a Banach space X, then G B{X) denotes the (bounded) inverse 
operator of T, with B{X) being the Banach space of bounded linear operators on 
X. Moreover, (j{T) denotes the spectrum of T, and T\y denotes the restriction of 
T to a subspace y of X. If Xi and X2 are two subspaces of X, then X1+X2 denotes 
their direct (but not necessarily orthogonal) sum. Bp{Ti) denotes the Schatten-von 
Neumann ideals of compact operators on a Hilbert space TH. with singular values in 

P > 1- A generic constant is denoted by c or C, and we also use the abbreviations 
M+ = [0,00), IR_ = (-00,0]. 



2. Preliminaries 

Throughout this paper we make the following assumptions on A and ii in (1.1) 
and (1.2): 
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Hypothesis 2.1. Suppose 

A,ReLl,(Rf-''. (2.1) 

Let $ denote the fundamental matrix solution of the unperturbed equation (1.1) 
so that ^'{x) = A{x)^{x), x gM. and $(0) = Id- Throughout, we assume that the 
propagator ^{x)^{x')~^ of (1.1) is exponentially bounded on M, that is, we suppose 
that 

sup ||$(x)$(a;')"lc<ix<i < oo, (2.2) 

x,x'eM, |x-a;'|<l 

or equivalently, there exist constants C G [1, oc) and a G R, such that 

||$(a;)$(a;')"^||c<i>«i < Ce"l^-^'l, a;,a;'GK. (2.3) 

First, we recall basic definitions and some standard facts regarding Bohl and Lya- 
punov exponents and dichotomies (see, e.g., [8, Scc.l], [10, Ch. III]). 
A projection Q in C'' is said to be uniformly conjugated by $ on M if 

sup||4>(a;)Q$(x)-^||cdxd < oo. (2.4) 
xeM 

The geometrical meaning of this condition is that the angle between the range and 
the null space of the projection ^{x)Q^{x)~^ in is uniformly separated from 
zero for all a; G M (see, e.g., [10, Corollary IV. 1.1]). 

Let Q be a projection in which is uniformly conjugated by $ on R. The 
upper Bohl exponent on R associated with the projection Q, denoted by >c[Q), is 
defined as the infimum of all x G R, such that for some constant C(x) G [l,oo), 
the following inequality holds for all a;, a;' G R such that x>x', 

\\^{x)Q^x'y\\cdxa < C(x)e^(^-^'). (2.5) 

The lower Bohl exponent on R, denoted by >c'{Q), is defined as the supremum 

of all >^ G R, such that (2.5) holds for all x < x\ where x,x' G R. 

Similarly, given a projection Q which is uniformly conjugated by $ on R_(_ (resp., 
on R_), the upper Bohl exponent, h+{Q), and the lower Bohl exponent, >c'^{Q), on 
R_i_ associated with the projection Q (resp., the Bohl exponents >c^{Q) and >c'_{Q) 
on R_), are defined in the same way except that one takes x, x' G R-|- (respectively, 
x,x' G R_) in (2.5) . The upper and lower Bohl exponents on R_|_ can also be 
expressed as follows (see [10, Theorem III.4.4]) 

fn\ V log\\^{x)Q<^>{x')-'^\\cdxd 
x+{Q) = limsup ; , (2.6) 

(x— x')— >oo, x'— >oo 

/ V log\\<i>{x')Q<^{x)-^\\cdxa 
x+{Q) = - limsup ; (2.7) 

{x — x') — ^00,2;' — *oo ^ 

similar formulas hold for x-{Q) and >r'_(Q). 

We introduce the upper, A+(Q), and the lower, A'_^(Q), Lyapunov exponents on 
R+ associated with the projection Q by the formulas 

\ fn\ V log||$(x)0||c<ixci , , log||Q$(a;)-i||c<ix<i 
A+ (Q) = hm sup — , A_,_ (Q) = — lim sup — — , 

X— »oo X X— ^00 X 

(2.8) 
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and the Lyapunov exponents on M_ by the formulas 

\ tr^\ V log||g$(a;)-i||c^xd , log||$(x)Q||c<ix 
A_((5)=limsup J— ^ — , A_ (Q) = - lim sup ■ 



\x\ x—* — oo \x\ 

(2.9) 

We remark that the upper (lower) Bohl exponent measures the best (worst) 
exponential growth of the propagator $(x)$(x')~^ and the upper (lower) Lyapunov 
exponent measures the best (worst) exponential growth of the fundamental solution 
$ relative to the projection Q. Under the assumption that the propagator of (1.1) 
is exponentially bounded, the Bohl exponents are finite [10, Theorem III. 4. 2]. We 
note the inequalities 

k'{Q) < k'^{Q) < A'±(Q) < A±(Q) < K±{Q) < k{Q), (2.10) 

and stress that they are strict in general (in particular, the possible inequality of 
the Bohl and Lyapunov exponents is exhibited in the classical Perron example, see, 
e.g., [10, Sect. III.4.4]). 

Remark 2.2. The set of Bohl and Lyapunov exponents for the resettled equation, 

y'{x) = (A{x)-nId)y{x), xGR,fieC, (2.11) 

with the associated propagator e~'^(^~^ )$(x)$(a;')~^, x,x' S R, is obtained from 
the set of the exponents of equation (1.1) by shifting the latter by the amount 
-Re(/i). O 

Example 2.3. If the unperturbed equation (1.1) is autonomous then the Bohl and 
Lyapunov exponents are the same and equal to the real parts of the eigenvalues of 
A. To give more details and fix notations, suppose that A{x) = A, x gW, for some 
A G <C^^'^, and let G C, fc = 1, . . . , d, be the eigenvalues of A. We split the set of 
the eigenvalues as (j{A) = U^^^Sj, where 1 <d' <d and for each j all eigenvalues 
that belong to have the same real part denoted by Hj, j = 1,. . . ,d'. Let Qj 
denote the spectral projection for A such that o'(j4|ran(Q )) = '^j- Then the Bohl 
and the Lyapunov exponents on M. (and on IR.±) associated with Qj, are equal to 
Xj (we note that we will always enumerate in increasing order of magnitude, 
X,- < Hj' for j < f, j, f = l,...,d'). O 

Returning to the general nonautonomous case, equation (1.1) is said to have 
a bounded dichotomy Q on IR. if Q is a projection in (called the dichotomy 
projection) so that for some constant M e [l,oo), the following inequalities hold 
for all x,x' eR: 

||$(.t)Q$(x')^'||c-x. <M, x>x', 

||$(.x)(/d-Q)$(a;')"'llc^x. <M, x<x'. 

Equation (1.1) is said to have an exponential dichotomy Q on M if Q is a projection in 
so that for some positive constants x and x', and some constants C(x), C(x') e 
[1, oo), the following inequalities hold for all x, x' e M: 

||$(a;)Q$(x')-'llc^x. < C(x)e-"("-"'\ x>x', (2.13) 

||$(x)(/d-Q)$(x')-ic^x. < C(x')e-'(---'), x<x'. (2.14) 

We note that Q is an exponential dichotomy for (1.1) on R if and only if the 
following inequalities for the Bohl exponents hold: 

x(Q) < < x'(/d - Q). (2.15) 



EVANS FUNCTIONS, JOST FUNCTIONS, AND FREDHOLM DETERMINANTS 



9 



Dichotomies on R+, respectively on M_, are defined in the same way except in 
(2.12) or (2.13), (2.14) one takes only x,x' G M+, respectively, x,x' G M_. 

Example 2.4. Assume in Example 2.3 that a{A) n = and d' > 2. Consider 

ko G {1, . . . , rf' — 1} such that >Ck„ < < Xko+i- Then Q = J^'j'LiQj is the 
exponential dichotomy projection for the equation y'{x) = Ay{x) on R, R_|_, and 
M_. O 

Definition 2.5. ([10, Sect. IV.4].) Assume d > 2, and let 1 < d' < d. A system 

{Qj}j=i of disjoint projections in is called an exponential splitting for (1.1) on 
R if the following conditions hold: 

(i) Ei=i Qj = Id- 

(a) Each projection Qj is uniformly conjugated by $ on R. 
{Hi) The segments [x'((5j), x(Qj)] are disjoint. 

The segments [x'{Qj), x{Qj)] of the real axis arc called the Bohl segw.ents as- 
sociated with Qj] they are determined by the lower and upper Bohl exponents 
on R associated with the projections Qj, and in what follows they will always be 
numbered so that >c{Qj) < >c'{Qj+i), j = 1, . . . , d' — 1. If d' > 2 and (1.1) has an 
exponential dichotomy Q on R, then for some fco, 1 < fco < — 1, we have the 
following splitting: 

ko d' 

Q = Y,Qj &nd h-Q= (2-16) 

3 = 1 3=ko + l 

SO that >c{Qka) < < V(Qfeo+i). 

Replacing R by R+ or R_ in Definition 2.5, one can also consider exponential 
splittings on R_|_ or R_ . We note that the exponential dichotomy projection Q on 
R is uniquely defined, while exponential dichotomy projections on R_|_ or R_ are 
not. Indeed, for example, in the case of M+, the subspace of initial data of 
bounded solutions of (1.1) on R+ is determined uniquely; to obtain an exponential 
dichotomy projection Q on R_|_ one can choose an arbitrary direct complement of 
the subspace Ng in so that = Ng+N^, and define the exponential dichotomy 
projection Q such that ra,n{Q) = Ng and ker(<5) = 7V„ (see [10, Remark IV.3.4]). 
We recall that if Q is an exponential dichotomy for (1.1) on R, then Q is also an 
exponential dichotomy on both semi-axes R+ and R_. Similarly, an exponential 
splitting {Qj}'j^i on R induces exponential splittings on R+ and R_; we stress that 
in view of (2.10), the Bohl segments for R are, generally, wider than the segments 
for R_|_ and R_ . This may happen since an exponential dichotomy on R_|_ does not 
necessarily imply an exponential dichotomy on R. 

Example 2.6. In Examples (2.3) and (2.4) each Bohl segment [}<' (Q j) , >i(Q j)] 
associated with Qj degenerates into a single point >Cj, j = 1, . . . ,d' , while the Bohl 
segments associated with the (exponential dichotomy) splitting {Q,Ici — Q} are 
given by [Ki,>Ckg] and [Kk„+i, Xd']; we recall that in this case 

XI = inf{Re(i/) G R | J/ G a(^)}, xu^ = sup{Re(j/) G R | G (j{A), Re{v) < 0}, 

Xko+i = inf{Re(i/) G R | G cr{A), Re{u) > 0}, Xd' = sup{Re(z/) G R | G cr(A)}. 

(2.17) 

O 
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Remark 2.7. The Lyapunov exponents for nonzero C'-valued solutions y — y{x) 
of (1.1) are defined as follows: 

A±(y)=limsupi^^CM*l, A'±(y) = liminfi^^I^M*). (2.18) 

For any nonzero C*- valued solution y of (1.1), and for any given exponential split- 
ting {Qk}t=i^ there exists a j € {1, ... , d!) such that all four Lyapunov exponents 
\±{y), )^'±{y) belong to the jth Bohl segment. One verifies that the actual jth. Bohl 
segment depends on the initial data yo = 2/(0) ^ of the solution y of (1.1), say, 
on 1R_|_, as follows: 

If yo e ran f ^ Ofc j then \+{y) < X+iQj), 

if t/o e ran ( ^ Qfc j \ ran ( ^ Qi) then A+(y) > X+{Qj), 

if A+(y) < X+{Qj) then yo e ran QfcY (2.19) 

if X'+{y) > X+{Qj) then y{0) G ^Xranfj^Q')^ 

where j G {1, . . . , d'} and we set Qo = 0. Similar assertions hold for the Laypunov 

exponents on ]R_ . O 

For future references, we record the assumptions on the coefficient A of the 
unperturbed equation (1.1). 

Hypothesis 2.8. Assume that A G Ll^^{Ry^'^. In addition, we suppose that the 
propagator ^{x)^{x')~^ , x,x' € M, is exponentially bounded on M., that is, 

sup \\^{x)^{x')~^\\cdxd < oo, (2.20) 

XjX'eK, \x-x'\<l 

and that (1.1) has an exponential dichotomy Q onM.. 

Turning to the perturbed equation (1.2), we assume that the perturbation R in 
(1.2) satisfies the condition 

\\R\yxaeL\R). (2.21) 
Let U and \R\ denote the dx d matrices in the polar decomposition of R: 

R{x) = U{x)\R{x)\, \R{x)\ = {R{x)*R{x)y/'^, xgM. (2.22) 
Throughout this paper we will use the notation 

Ri{x) = U{x)\R{x)\i, Rr{x) = \R{x)\i , x gR, (2.23) 

so that 

R{x) = Re{x)Rr{x), xeR. (2.24) 
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Lemma 2.9. Assume Hypothesis 2.8 and condition (2.21). Then the integral op- 
erator K with integral kernel given by 

''^'''''>-\Rr{xMx){I,-QMxr'Rdx'), x<x', ^'-"'^ 
^5 a Hilbert-Schmidt operator on L^{M.Y. Moreover, the map 



ii(]R)rfxd^g2(L2(M)rf) 



(2.26) 



^5 continuous. 



Proof. Since by hypothesis equation (1.1) has an exponential dichotomy Q on R, 
the projection Q is also a bounded dichotomy on R, and thus (2.12) holds. Using 
(2.12), the Hilbert-Schmidt norm of the integral operator K can be estimated as 
follows: 

\\K\\B2{L^Ry)= / dxdx' \\K{x,x')\\ldxd 

= [ dx f dx'\\Rr{x)<l>{x)Q<^{x')-^Re{x')\\ld.4 
+ [ dx f dx'\\Rr{x)^x){Id - Q)^x')-^Re{x')\\laxa 

JR Jx 

'■ dx da;'||i?(a;)||c<ixd||i?(a;')||c<ix<i 
Jr Jr 



< M 



= M yj dx\\R{x)\\cdxdj . (2.27) 
In the last estimate wc used the inequalities 

||i?r(x)||c^x. < ||i?(,x)||J,,„ \\RfXx)\\c^.. < \\R{x)\\l,,,, X&R, (2.28) 

which follow from the spectral theorem for self-adjoint matrices. The continuity 
of the mapping in (2.26) holds because the mappings R Rr and R i-^ Re are 
continuous from the Banach space Z/^(M)''^'' to the Hilbert space i^(]R)'^^''. □ 

Remeirk 2.10. Since the propagator of (1.1) is exponentially bounded, the formula 

{Elu){x) = ^{x-t)u{x-t), ueL^iRf, xeR,t>0, (2.29) 

defines a strongly continuous semigroup, {i?^}t>o, on I/^(M)'^, called the evolution 
semigroup (sec [7, Sect. 3.2] for a detailed discussion and extensive bibliography). It 
is well-known (see, e.g., [7, Theorem 3.17]) that the generator, G, of the semigroup 
in (2.29) has a bounded inverse if and only if (1.1) has an exponential dichotomy 
Q on R and, if this is the case, the inverse operator is an integral operator on 
L^(R)'' with integral kernel given by 

^^^'^^-\$(.)(/,-0)$(x')-\ x<x', ^'-'"^ 

where Q is the dichotomy projection. Thus, under the hypotheses of Lemma 2.9, 
the integral operator K admits the representation K = Mr^G~^M.r^, where Ain^ 
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and denote the maximally defined operators of multiplication by the matrix- 

valued functions Rr,R£ G L^(R)'^^'', and the operator S = Mb.^G^^Mb., (with 
the bounded closure S) is defined on the set of functions u E dom{M r, ) such that 
e L2(K)dxd_ -y^g [27] for a detailed discussion of Bir man Schwinger-type 

operators and the Birman-Schwinger principle. O 

Remark 2.11. Using the semigroup in (2.29) one can describe the exponential 

dichotomy projection Q for (1.1) on R by means of the Riesz spectral projection 
for the operator E\,, with t = 1 in (2.29), as follows (see [7, Theorems 3.13, 3.17]): 
Assume that the generator G of the semigroup in (2.29) is invertible on i^(R)'^. 
By [7, Theorem 3.13], it follows that E]^ has no spectrum on the unit circle. Let 
Q denote the spectral projection for E\ on L^(R)'' such that cr(-E^$ |ran(Q)) is the 
part of the spectrum of E\, contained inside the unit disc. Theorem 3.17 in [7] 
states that Q is an operator of multiplication in i^(R)'^ by a projection- valued 
function Q{-) G Cb(R)''^'' and, in addition, Q = (3(0) is the exponential dichotomy 
projection for (1.1) on R. O 

Remark 2.12. The proof of Lemma 2.9 shows that the assumption of an exponen- 
tial dichotomy on R for (1.1) can be relaxed to require only a hounded dichotomy. 
Under this weaker assumption the generator G of the semigroup in (2.29) may 
not be invertible (cf. Remark 2.10), or may not even be a Predholm operator on 
L^(R)''. Indeed, by a well-known Dichotomy Theorem (sometimes called Palmer's 
theorem), the operator G is Fredholm on L'^{KY if ^-i^d O'^iy if (1-1) has exponential 
dichotomies Q- on R_ and Q+ on R+; see [3], [48], [49], [56] or [57, Theorem 3.2], 
and also [41], [53] for more recent versions of the dichotomy theorem. O 

Next, we will discuss the Bohl exponents and exponential splittings for the per- 
turbed equation (1.2). Since the propagator of (1.1) is exponentially bounded, 
and condition (2.21) holds, it follows from the variation of constants formula and 
Gronwall's inequality that the propagator of the perturbed equation (1.2) is also 
exponentially bounded (see, e.g., [33, Lemma IV. 4.1]). Thus, the Bohl exponents 
for the perturbed equation (1.2) are finite. 

We conclude this section with two technical results to the effect that first, 
the exponential dichotomy and exponential splitting of (1.1) on R+ persist un- 
der i^(R+)''^''-perturbations, and second, that the corresponding Bohl exponents 
do not change under these perturbations (similar facts hold for R_). A proof of 
the first part of these results can be found in [8, Proposition 4.1] or [10, Theorem 
IV.5.1]. However, we were not able to find the second part in the literature, and so 
we will briefly sketch a proof in Appendix B by modifying some of the arguments 
in [10, Ch. IV]. Lemma 2.13 is used in Section 6, while its direct generalization. 
Lemma 2.14, is used in Section 7. 

Lemma 2.13. Assume that the propagator of the unperturbed equation (1.1) on 
R+ is exponentially bounded, that the unperturbed equation (1.1) has an exponential 
dichotomy Q on R+, and the perturbation R satisfies ||i?||c<ix<! e L^(R+). Then the 
perturbed equation (1.2) has an exponential dichotomy P on R+ such that 

x+{P) = x+{Q) and >t'+{h - P) = x'+ih - Q), (2.31) 

where >c+{P) '^'^'^ ^^'+{1(1 — P) (I're the Bohl exponents for the perturbed equation 
(1.2) on R+. Here, P is the projection in parallel to ker{Q) onto the subspace 
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consisting of the values y{Q) at zero of all bounded solutions y of (1.2) on M+. 
In addition, 

dim(ran(P)) = dini(ran(Q)) and dini(ker(P)) = dini(ker(Q)). (2.32) 

Lemma 2.14. Assume that the propagator of the unperturbed equation (1.1) on 
is exponentially bounded, the unperturbed equation (1.1) has an exponential splitting 
{Qj}j=ij ^ < d' < d, on M+, and the perturbation R satisfies ||i?||cdxd G L^(M_|_). 
Then, the perturbed equation (1.2) has an exponential splitting {Pj}j^i on IR+ such 
that 

dim(ran(Pj)) = dini(ran((5j)), j = 1, . . . ,d' , (2.33) 

and 

K+{P,) = K+{Q,), j = l,...,d', (2.34) 
where x+{Pj) are the Bohl exponents for the perturbed equation (1.2) on M+. 

3. Matrix- Valued Jost Solutions 

In this section we discuss optimal conditions on the perturbation R under which, 
assuming an exponential dichotomy of (1.1), one can establish existence and unique- 
ness of bounded matrix-valued solutions of (1.2) (we will call them matrix-valued 
Jost solutions) that are asymptotically close with respect to an exponential weight 
factor to bounded solutions of the unperturbed equation (1.1). These results are 
motivated by the study of the one-dimensional Schrodinger equation on the real 
axis, 

-u"{x) + V{x)u{x) = k^u{x), xeR, fc G C, (3.1) 
with an integrable potential V e L^{R;dx), where the Jost solutions, u±{k, •), are 
introduced as the solutions of (3.1) that arc asymptotically close to the free plane 
waves: u±{k,x) ~ g±ikx ^g^^^ ^ g^^ jg^ XVII]). 

X — *±oo 

Definition 3.1. Under Hypothesis 2.8, d x d matrix-valued solutions Y± of the 

differential equation 

Y±{x) = {A{x) + R{x))Y±{x), X e M±, (3.2) 

with a locally integrable d x d matrix-valued function R, are called matrix-valued 
Jost solutions if 

hmsup logr+(^)-^(^)Qllc^xa ^ ^^^^^ ^3 3^ 

and 

j.^.^^ logr-(.)-^(.)(/,-Q)||cax. ^ _ (3_^) 

x— >— oo X 

To motivate this definition, we first note that the Lyapunov exponents X+{Q) 

and A'+((5) of the bounded matrix-valued solution $ Q of (1.1) on belong to the 
Bohl segment [k'_^{Q), k_^.{Q)]. If F_|_ is a matrix-valued solution of (1.2) satisfying 
(3.3), then 1+ is bounded and therefore Lemma 2.13 implies that its Lyapunov 
exponents belong to the same segment. Thus, the significance of (3.3) is that the 
solution approximates the solution $ Q exponentially better than either of these 
solutions decays as a; — > oo. 

Our first result shows the existence and uniqueness of the matrix-valued Jost 
solution 1+ on using a rather strong assumption on the exponential fall-off of 
the perturbation R. 
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Theorem 3.2. Assume Hypothesis 2.8 and the condition 

|li?llccixd e Li(R+;e'^^(ia;) (3.5) 

for some constant /3 > satisfying the inequality 

a > X+{Q) - x'4Q). (3.6) 

Then the perturbed differential equation (1.2) has a unique matrix-valued Jost so- 
lution Y+{x), X > 0. Moreover, the map 



i?t->y+(0) ^^'^^ 



is continuous. 



Proof. We split the proof of Theorem 3.2 into three steps. 

Step 1. Uniqueness. Suppose that Y+^i and 1+^2 are two different dxd matrix- 
valued Jost solutions of (1.2), and introduce the nonzero solution Y(x) = y+,i(a;) — 
y+,2(a;), a; > 0. Using (3.3), we infer: 

log||r(a;)||c<ixd / log||y+,fe(a;) - $(a;)Q||cdxd\ 
limsup — < max I limsup I < Hj^{Q). 

X— ^oo ^ fc=l,2 Y X— »oo ^ J 

(3.8) 

This shows, in particular, that F is a bounded solution of (1.2). Since ||i?||c<ix<i € 
i^(K+) by hypothesis, Lemma 2.13 yields an exponential dichotomy P for (1.2) 
such that >i'j^{P) = x'+iQ). Since F is a bounded solution, F(0) € ran(P) and 
hence we arrive at the contradiction 

hmsupi^^rak^ll > V+(P) = x;(g), (3.9) 

X— *oo X 

proving the uniqueness part of the theorem. 

Step 2. Existence. By hypothesis, (1.1) has an exponential dichotomy Q on M_|_. 
Therefore, for any e > there exists a constant C{e) >1 such that applying (2.5) 
in the definitions of x'j^{Q) and >c'^{Id — Q), we obtain the estimate 

mxMx'r'Wc... < mx)Q^xT'\\c^x^ + mx){id - qmx')-'\\c^>'^ 

where we denoted 

f,{x) = C(e)e^^ a; > 0. (3.11) 
Since x^(Q) < x'_^{Id — Q), we get the estimate 

\\^x)^{x')-^\\cdxi < 2e^+('5)(^-^'Vc(k - x'\), {)<x<x'. (3.12) 
Next, fix £ such that < 2£ < /3 — (A+(Q) — x'_^{Q)), introduce the function 

p(a;) = e(^+('3)--V(Q))-/|(a;)||ii(a;)||cdx., ar > 0, (3.13) 

and observe that < p G i^(K+) due to (3.5) and (3.6). 

On the space Cb(IR+)'^^'' we now define a Volterra integral operator, V, by 

POO 

{VZ){x)= / dx' e-^+^'^^^''-'''^fe{x)-^fe{x')^{x)^{x')-^R{x')Z{x'), (3.14) 

J X 

and consider the corresponding Volterra integral equation 

Z{x) = Z(^\x)-{VZ){x), x>Q, (3.15) 
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where Z^^^ is defined by 

Z(o)(a;) = e-^+(«(")/,(x)~i$(x)0, a; > 0. (3.16) 

Using (3.12) and the inequality >c'_^{Q) — ^+{Q) < (cf. (2.10)), one derives the 
following estimate for the integral kernel of the integral operator V: 

||y(x,ar')||c^x^ < 2eK(Q)-A+(Q))(--')/^(a;)-i/^(a;')/^(|a; 

<2eKW)-^+(Q))-/,(a;)-ie-(-+(«)-^+(«))^'/2(^/)||^(^/)||^^^^ 

= 2e(^+(«)-^+(«))^/e(a;)-V(a;') (3.17) 
< 2fs{x)-'^p{x'), 0<x<x'. (3.18) 

Since fe{x) > 1, x G R+, and p e L^(M-i-), the estimate (3.18) shows that F is a 
compact operator on Cb(IR+)'^^'' with spectral radius equal to zero. Since the expo- 
nential decay of H^QUc^xd is controlled by X+{Q), we see that G Cb(]R+)''^''. 
Thus, the integral equation (3.15) has a unique solution Z G Cb(K+)'^^'^ that can 
be obtained by the iteration process 

oo 

Z{x) = Z'^^Hx), Z'^\x) = (yz(^'-i))(a;), j G N, a; > 0. (3.19) 

Moreover, from (3.15) and (3.17) it also follows that 
Mx)\\Z{x) - Z^''\x)\\c<i.^ 

POO 

< 2||Z||c,(M,)-.eKW)-^+W))- / dx'pix'). 

A straightforward computation shows that the function 

Y+{x) = e/+(«-/,(x)Z(x), x > 0, (3.21) 

yields a solution of the perturbed differential equation (1.2), and hence (3.20) and 
p G L^(R+) imply 

e-<iQ>\\Y+{x)-^{x)Q\\c^.^ = o(l). (3.22) 



(3.20) 



X — ^oo 



In addition, using (3.5) and (3.6) again, we may choose (5 > sufficiently small 
such that the function P5{x) = e''-'^p{x) is integrable on R+. If a; < x' , then 
p{x') < e~^^ps{x'), and thus (3.20) implies 



Mx)\\Z{x)-Z(''\x)\\c^.. 

<2||Z||c,(M+).x.eKW)-^+(Q))-e-*- / dx'psix'). 

J X 



(3.23) 



This gives a better asymptotic relation than (3.22), 

e-<(Q)-\\Y+{x)-^x)Q\y.. = o(e-^^), (3.24) 

leading to 

log\\Y^{x) - Hx)Q\\e... ^ ^^^Q^ _ ^ ^ ^^(^^^ (3 25) 
proving the existence part of the theorem. 
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Step 3. Continuity. Let {-Rnj^i be a sequence of d x d matrix-valued functions 
converging to R in L^(R+; e'^^dx)''^'^ as n — > oo. On the Banach space Ch{^+Y^'^ 
introduce Voltcrra integral operators Vn by 

/•oo 

{VnZ){x)= / dx'e-^+('^^^''-'''^fe{x)-^h{x')^{x)^{x')-^Rn{x')Z{x'), n € N. 

J X 

(3.26) 

An estimate similar to (3.18) shows that the sequence of the operators Vn converges 
in operator norm to the integral operator V, which in turn, yields convergence in 
Cb(]R+)''^'^ of the unique solutions Z„ of the Volterra integral equations 

Z^{x) = Z^'^^x) - {VnZn){x), X>Q, (3.27) 

to the unique solution Z of the integral equation (3.15). In particular, 

lim Z„(0) = Z(0) in C^""^ and hence lim F+,„(0) = F+(0), (3.28) 

where, similarly to (3.21), we denote F_|__„(a;) = e^+^^^'' fe{x)Zn(x). This proves 
continuity of the mapping (3.7). □ 

Remark 3.3. This theorem holds under a weaker assumption than the exponential 
dichotomy hypothesis (cf.(2.15)). In fact, it is sufficient to require that x+((3) < 
><M-(-^d - Q) and x+((5) < only. O 

Remark 3.4. Replacing assumptions (3.5) and (3.6) by the hypothesis 

||ii||c<ix<i e L\R+;e^''dx) with /3 > -X'+{Q), (3.29) 

the matrix-valued Jost solution can be uniquely determined by solving the 
Volterra integral equation 

/•oo 

Y+{x) = ^{x)Q- dx' ^{x)^x')-^R{x')Y+{x'), x>Q; (3.30) 

J X 

that is, in this case, the "modified" Volterra integral equation (3.14), (3.15) on 
Chi^+Y^"^ i'' i^ot needed. Indeed, assume (3 > —X'_^_{Q), choose e such that < 
£ < min{/3 + X'_^_{Q) , — Q)}, and use the inequalities x' > a; > and X'+{Q) < 
< — Q) to estimate the integral kernel of the integral operator in (3.30) as 
follows: 

\mx)^x')-^R{x')\\c^., < {\mx)Q\\c^.4Q^x')-^\\c^., 

+ mx){I - Q)$(x')-'||c<^x.) \\R{x')\\c^.. 

' (3.31) 

< (C(e)e(-^V(0)+^)^' +C(e)e(^+('f-'3)-^)(^--'))||ii-(x')||c^x. 

< 2C7(e)e-(^+(«+'^^-")"'e'^^"-'||i?(,x')||c^x. < 2C{e)e<^^' \\R[x')\\cd-.d. 

We recall that x^(Q) < AY((3) < A+(Q) < 0. Condition (3.6) is formulated 
in terms of the difference A+(Q) — >c'^{Q) while condition (3.29) is formulated in 
terms of the quantity —X'_^_{Q). These two numbers, A+(Q) — >c'^{Q) and — A^((5) 
are fairly independent, and thus either of the two conditions, (3.6) and (3.29), can 
be better than the other. Thus, the existence of matrix- valued Jost solutions follows 
under the assumption /? > min{A+((3) — x^((5), —X'_^_{Q)}. O 
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Remark 3.5. Example 3.6 below shows that the statement of Theorem 3.2 is opti- 
mal in the sense that, under the exponential dichotomy hypothesis, the assumption 
in (3.6) on the exponential fall-off of the perturbation R in general cannot be re- 
laxed: Indeed, if condition (3.6) is violated, then the matrix-valued Jost solution 
may not exist. Example 3.6 also shows that if Definition 3.1 is modified to replace 
the strict inequality in (3.3) by the condition 

hmsup loglin(^)-^WIIc^x^ < (3 32) 

then the matrix-valued Jost solution in this modified sense is generally not unique. 

O 

Example 3.6. Consider equations (1.1) and (1.2) with 

0\ /O XM+(a;)e-^cosa; 0\ 

A{x) = I -1 and R{x) = , a; G M, 

1/ \0 0/ 

(3.33) 

where xm+{x) denotes the characteristic function of = [0, oo). The fundamental 
solution <1>, the exponential dichotomy Q on M, and the bounded solution f&Q of 
(1.1) on IR+ are given as follows: 

0\ 


xgR. (3.34) 

The Bohl and Lyapunov exponents on R_)_ associated with the dichotomy projection 
Q are x'+iQ) = X'+{Q) = -2 and ><+(Q) = X+{Q) = -1. Clearly, \\R\y^^ € 
L^(IR+; e^'^dx) with /3 = X+{Q) — x'^{Q) = 1, and hence condition (3.6) is violated. 
A direct computation shows that the perturbed equation (1.2) has two linearly 
independent bounded C^-valued solutions on M_|_ 

(e~^^\ /e~^'^sinx\ 
j , y^{x) ={ e-^ 1 , X > 0, (3.35) 

and thus any bounded C'-valued solution of (1.2) on M_|_ is a linear combination 
of these two. Hence, any bounded matrix-valued solution Y of (1.2) on M_|_ is 
necessarily of the form 

/g-2x e-2^(sinx-hC) 0\ 
Y{x) = e-^ , a; > 0, (3.36) 

\ 0/ 

where C G C is an arbitrary constant. Thus, it satisfies 

hmsup l°grW-^WQIIc^x- = _2 = ^;(Q). (3.37) 

X — ^oo X 

This observation shows that the dx d matrix-valued Jost solution as introduced in 
Definition 3.1 does not exist if (3.6) fails, while the solutions satisfying (3.32) are 
not unique. O 
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Remark 3.7. In the proof of Theorem 3.2 we suggested a method of introducing 
the Jost solution by Y+(x) = e''^Z(x), x > 0, where Z e Cb(K+)'*'"^ is a solution 
of the "/i-modified" Volterra integral equation 

/■oo 

Z(x) = e-''^$(a;)Q - / dx' e-''^''-''">^x)^x')-'^R{x')Z{x'), x>0, (3.38) 

for an appropriate choice of ^ G R. We stress that this method fails if the expo- 
nential fall-off hypothesis (3.6) is violated. Indeed, for the equations in Example 
3.6, if ^ < X+{Q) = —1, the term e~^^^{x)Q in the integral equation (3.38) is 
not a bounded function, while if /x > x+{Q) = —1, the Volterra integral opera- 
tor is unbounded, both in the Banach space Cb(IR+)'^^'' and in the Hilbert space 

Similar to the proof of Theorem 3.2 (replacing >c'_^{Q) by >c-{Id — Q) and X+{Q) 
by \'_{Id — Q)) one proves the following result for the negative half-line M_. 

Theorem 3.8. Assume Hypothesis 2.8 and the condition 

\\R\\cd^d G L^(M_;e'^l^lrfa;) for some /? > x_(/d - Q) - \'_{Id - Q)- (3.39) 

Then the perturbed differential equation in (1.2) has a unique matrix-valued Jost 
solution Y-(x), x <Q. Moreover, the map 

I i?H^Y_(0) ^ ■ ^ 

is continuous. 

4. The Perturbation Determinant 

In this section we employ results from [28], succinctly summarized in Appendix 
A, to express the infinite-dimensional (modified) Fredholm determinant of the inte- 
gral operator I + K, given by (2.25), via a finitc-dimcnsional determinant obtained 
by means of the matrix-valued Jost solutions introduced in Definition 3.1. This 
result is a generalization of the classical relation between a Fredholm determinant 
and the Jost function due to Jost and Pais [36] (sec, also [46, Sect. 12.1.2]). As 
we will see in Section 9, the finite-dimensional determinant is related to the Evans 
function associated with (1.1) and (1.2). We recall Lemma 2.9, and also the gen- 
eral fact (see, e.g., [30, Ch. XIII], [31, Sect. IV.2], [65], [66, Ch. 9]) that if K is any 
Hilbert-Schmidt operator, then the 2-modified Fredholm perturbation determinant 
det2 (/ + K) is given by the formulas 

det2(/ + /i:) = det((7 + /i:)e-^) = JJ ((1 + A)e-^), (4.1) 

counting algebraic multiplicites of the eigenvalues of K. Moreover, if ii' is a trace 
class operator, then 

det2{I + K)=e-'^''^^^det{I + K), where det{I + K)= (1 + A). (4.2) 

Throughout, we will use the following notation: 

poo rO 

e= dxtvc''{^x)Q^x)-'^R{x)) - dxtTcdmx){Id-Q)^x)-'^R{x)). 

Jo J -oo 

(4.3) 
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Theorem 4.1. Assume Hypothesis 2.8 and the condition 

||i?||c^xd eLi(R;e''l"lcix) (4.4) 

for some 

/?>max{A+(0) ->^;(0),x_(/d-0)-A'_(/<i-0)}. (4.5) 

Let K be the integral operator on L'^iM)'^ whose integral kernel is given by (2.25). 
Then the 2-modified perturbation determinant det2{I+K) admits the representation 

det2(/ + i^) = e^detc^(F+(0) + y_(0)), (4.6) 

where & is defined in (4.3), and Y± are the matrix-valued Jost solutions on M±, 
respectively. 

Proof. Assume temporarily that 

supp (i?) is compact. (4.7) 

Then ||i?||c<ixd e Li(R;e''l^l(ix) for any real (3. In particular, by Remark 3.4 and 
its obvious IR_-analog, the d x d matrix-valued Jost solutions Y± are the unique 
bounded solutions of the following Volterra integral equations on R±, 



f'OC 

Y+{x) = $(x)Q - / dx' <^>{x)<^{x')-^R{x')Y+{x'), x>0, 

J X 

y_(x) = ^x){Id -Q)+ dx' ^x)^x')-'^R{x')Y_{x'), x<0. 

J — OO 



(4.8) 



We note that i + (a;) = Y+{x)Q and Y-{x) = Y-{x){Id — Q) by uniqueness of the 
solutions. Moreover, by (4.7) we have Yj^{x) = ^{x)Q iov x > n and Y-{x) = 
^{x){Id — Q) for X < ~n for sufficiently large n G N. 

Treating matrices as operators on respective spaces, we introduce the following 
notations for a; e M, 

fi{x) = Rr{x)^{x)Q : ran(Q) ^ 

f2{x) = Rr{x)^{x){Ia - Q) : ker(g) ^ C^, 

gi{x) = Q^x)-^Re{x) : & ^ ran(g), 

92{x) = -{Id - QMx)-^Re{x) : & ^ ker(g). 

In addition, let di = dim(ran((5)) and d2 = dim(ker((5)), and denote 

H{x,x') = fi{x)gi{x') - f2{x)g2{x') = Rr{x)^{x)^{x')-^Re{x'), x,x' e R. 

We note that the functions fj and gj, j = 1,2, arc compactly supported on M due 
to (4.7), and thus assumption (A.l) holds. Hence, the results recorded in Appendix 
A are at our disposal. Next, we introduce the Volterra integral equations 

/•oo 

fi{x)=fi{x)- dx' H{x,x')h{x'), xeR, 

■'}x (4.10) 
h{x) = f2{x) + / dx' H{x,x')f2{x'), xeR, 

and let fi{x) : ran((5) C and f2{x) : ker((5) be the unique solutions of 

(4.10) that satisfy fj e L^^R^dxd,^ j = 1,2 (cf. Appendix A). The dxd matrix 
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U{x), X e M, defined in (A. 9), written as a block-operator with respect to the 
decomposition = ran(Q)-i-ker(Q), then reads as follows: 

UM Cdx'gi{x')h{x') j:^dx'9^{x')Mx') \ ..... 

^^^^ 1^ Cdx'g2ix')Mx') h^-r-^dx'g,{x')h{x'))- ^^'''^ 

Writing R{x) = Re{x)Rr{x), and multiplying (4.8) by Rr{x) from the left, we arrive 
at (4.10) and observe that 

/i(a;) = Rr{x)Y+{x)\^^niQ), x>0, 

f2{x) = Rr{x)Y_{x)\i,er{Q), X < 0. 

Setting a; = in (4.8), using (4.9) and (4.12), and writing the matrix y+ (0) + F_ (0) 
as a block-operator with respect to the decomposition = ran(Q)-t- ker((3), we 
conclude from (4.11) at a; = that ;7(0) = y+(0) -^^-(0). Since the integral kernel 
K{x, x') in (2.25) can be represented as 

y-T2{x)g2{x'), x<x', 

Theorem A.l (m) proves (4.6) for compactly supported perturbations. Specifically, 
one can use (A. 13) with a;o = and then apply some elementary properties of 
matrix traces. 

To remove the compact support assumption on the matrix-valued function R 
one proceeds as follows. Given n e N, we introduce the truncations 



R^^:,) = ry>^ (4.14) 




Let Rn{x) — Vn{x)\Rn{x)\, x G M, bc the polar decompositions of i?„(x), denote 
Re,n{x) = Vn{x)\Rn{x)\^l'^ and Rr^nix) = \Rn{x)\^l'^ ^ and introduce the integral 
kernels of Hilbert-Schmidt integral operators on by 



K (x x'^ - I -RrAxmx)QHxT'RiAx'). X>X', , ^ .4 ^5. 

' ' ~ ^ RrAx)Hx){h - Qmx'r'R^Ax'), x<x\ "^'"^ ^ ^- ^^-'^^ 



Since the support of i?„ is compact, one infers by the first step of the proof that 
det2(/ + i^„) = e®"detcd(r+,„(0) + Y_,„(0)), n e N, (4.16) 

where 



/■oo 

e„= / dxtrc<i{^{x)Q^{x)-^Rn{x)) 
Jo 

dxtTc4Hx){Id - Q)<P{x)-^Rnix)) 



(4.17) 



00 



and Y±^„ are the matrix- valued Jost solutions of the tnmcated perturbed equation 
y' = {A + Rn)y on Since i?„ converges to R in L^(M)''^'' as n ^ cx), one 

concludes that 

lim e„ = e, (4.18) 
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where Q is given by (4.3). Using the estimate in the proof of Lemma 2.9, one checks 
that the operators Kn converge to the operator K in Hilbert-Schmidt norm, and 
hence (see, e.g., [30, Ch. IX], [66, Ch. 9]) 

Hm det2(/ + K^) = det2(/ + K). (4.19) 

n— >oo 

Finally, since the sequence i?„ converges to R in L^(R; e'^'^'dx)''^'' as n — > oo, 
applying assertions (3.7) and (3.40) of Theorems 3.2 and 3.8, we obtain 

lim y±,„(0) =y±(0). (4.20) 

n — -'Oc 

Combining (4.18)-(4.20) then completes the proof. □ 

Remark 4.2. The following heuristic argument may be helpful in understanding 
the role of the factor e® in (4.6) and shows that the appearance of this factor is 
quite natural. First, we observe that the integral kernel (2.25) of the operator K is 
generally discontinuous on the diagonal x = x'. Accordingly, the definition of 8 in 
(4.3) can be re-written as follows: 

POO 

-Q= dx trcd{K{x + 0,x)) + dx ty:cd{K{x - 0,x)). (4.21) 

Jo J-oo 

Heuristically, the right-hand side of (4.21) can be viewed as a "regularized integral 
trace" of the operator K. Next, we assume in addition that K is a trace class 
operator (this requires quite a stretch of imagination!), and that tr(i^) equals the 
integral trace — O, 

tr(if ) = -e. (4.22) 
Then formally applying (4.2) and (4.6) yields 

det2(/ + K) = e-*''(^) det(/ + K) = e®detc<j(y+(0) + r_(0)) (4.23) 
and hence formally, 

det(/ + K)= dctcd(r+(0) + r-(0)). (4.24) 

As we will see in Section 9, the determinant detc<i(i + (0) + Y-{0)) is in fact the 
Evans function for the equations (1.1) and (1.2). <> 

Remeirk 4.3. The proof of Theorem 4.1 shows that the 2-modified determinant of 
the operator I + K can be computed by the formula 

det2(/ + K) = e® lim detc<i(F+,„(0) + F-,„(0)) (4.25) 



n—>-oo 



under the weaker assumption ||i?||cdxd G L^{R) as opposed to (4.4) and (4.5). Here 
y±,n are the unique matrix-valued Jost solutions of the Volterra integral equations 

POO 

Y+,„{x) = ^x)Q - / dx' ^x)^x')-^Rn{x')Y+,n{x'), x>0, 

p. (4.26) 
F_,„(x) = $(x)(/d - Q) + / dx' $(x)$(x')~'i?„(x')r_,„(x'), x<0, 



and Rn is the truncated perturbation (4.14). We also have the following identities: 

F+,„(o) = y+,„(o)Q, y_,„(o) = y-,„(o)(7d - q), (4.27) 

since the solutions Y± „ are unique. O 
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Remark 4.4. The exponential fall-off assumption (4.5) has been imposed in Theo- 
rem 4.1 to make sure that the matrix- valued Jost solutions Y± exist and are unique 
(cf. Remark 3.5). Moreover, under this assumption, one uses continuity of the 
mapping (3.7) to conclude that 



Y. 



■,n(0))) 



(4.28) 



lim detc<i(y+,„(0) + y-,„(0)) = detc. lim (r+,„(0) + 

n — >-oo \ n — >-oo 

= detc.(y+(o) + r_(o)). 

Example 3.6 shows that the exponential fall-off assumptions (4.5) are indeed needed 
to interchange the limit n — > oo and the determinant in (4.28). In fact, in this 
example one immediately verifies that the Jost solutions 1±,„ associated with the 
truncated perturbation i?„ arc given by 





g-2. 


sin 


{x)- 











V 























•) 




V 























a; < 0. 












sm n 



)]) 



< X < n. 



X > n, 



(4.29) 



Clearly, the sequence of the 3x3 matrices y+,„(0) does not converge as n ^ oo, 
while lim„^(x) detc'i(5:+,n(0)-|-Yl,n(0)) does exist since detc<i(^+,ra(0)-fYL,„(0)) ~ 1 
for all n G N. Consequently, (4.28) fails. O 

5. Sub-exponential Weights and the Jost Function 

As we have seen in Section 3, the exponential fall-off hypotheses (3.6) and (3.39) 

on the perturbation R cannot be relaxed in general because the unperturbed equa- 
tion might have at least one of the following two properties: First, the upper 
Lyapunov exponent associated with the dichotomy projection Q may not coincide 
with the lower Bohl exponent and then, necessarily, the corresponding Bohl interval 
is of positive width; and, second, the estimate (2.5) used in the definition of the 
lower Bohl exponent x'(Q) may not hold for x = >^{Q) but only for k = x'{Q) — e 
for any e > 0. This results in the de facto presence of an exponential weight factor 
fe{\x — x'\) = C{e)e^^^~'^ I on the right-hand side of (2.5); cf. also (3.11) and the 
effect of this on the proof of Theorem 3.2. In this section we restrict our attention 
to the class of unperturbed equations (1.1) that do not have either of these proper- 
ties. That is, we consider the case when each of the Bohl segments associated with 
the dichotomy projections Q and (7^ — Q) degenerates into a single point and, in 
addition, we assume that the exponential weight factors in the estimates such as 
(2.5) are replaced by sub- exponential weight factors induced by a given monotone 
weight function /. 

We introduce a weight function, f, which, by definition, is a nondecreasing func- 
tion / : M_|_ M_|_, satisfying the following conditions: 

log(/(a;)) 



lim sup ■ 



= 0, /(0)>1. 



(5.1) 
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Hypothesis 5.1. Assume Hypothesis 2.8. Suppose that the Bohl segments on M 
associated with the projections Q and {Id — Q) have zero width, that is, 

x'iQ) = x{Q) < < x'(7d - Q) = x{Id - Q). (5.2) 

Given a weight function f, we asssume, in addition, that the following estimates 
hold for all x,x' G R; 

||$(x)Q$(x')"^||c^xd < e^(^)(=""^'\/(|a; - x'\), < x < x' , 

Mx){Id - QMx')-^\\c<ix^ < e"(^^-«)("-"')/(|x - a;'|), > a; > a;'. ^^'^^ 

Given Hypothesis 5.1 we now introduce th.e matrix- valued Jost solutions as solu- 
tions of certain "modified" Volterra integral equations (cf. (3.14)-(3.15) and (4.8)). 

Definition 5.2. Assume Hypothesis 5.1 and suppose 

\\R\\c^.,eL\m;f^{\x\)dx). (5.4) 

Then matrix- valued solutions Y± of (1.2) on R± are called matrix-valued Jost so- 
lutions of (1.2) if 

F+(a;) = /(M)e-('3)-Z+(x), a: > 0, 

y_(a;) = /(|a;|)e"(^'^-'3)^Z_(a;), a; < 0, 

where Z+{x), x > 0, and Z-{x), x < 0, are the unique bounded solutions of the 
Volterra integral equations 

Z+{x) = e-''('5)^/(a;)-i$(a;)Q 

/>00 

- dx'e-''^'^'>^'=-''">f{x)-^f{x')<P{x)<P{x')-^R{x')Z+{x'), x>Q, 

J X 

Z_{x) = e-^(^<*-'^)^/(|a;|)-i$(a;)(7 - Q) (5.6) 
+ r dx' e-''^''-^^^'^-'''^ f{\x\)-^ f{\x'\)^{x)^{x')-^R{x')Z-{x'), x < 0, 

J — OC 

on M+ and M_. 

Almost literally repeating the proof of Theorems 3.2 and 3.8 (replacing /g by /), 
one concludes that the matrix-valued Jost solutions Y± in the sense of Definition 
5.2 satisfy the following asymptotic relations (cf. (3.22)): 

e--iQ>\\Y+{x)-^x)Q\\c^.. = 0(1), 
^->^(ia-Q)^Y_{x)-^{x){h-Q)\\c^.. = o{l). 
Moreover, one obtains the following version of Theorem 4.1. 
Theorem 5.3. Assume Hypothesis 5.1, suppose that 

\\R\\c^..eL\^-f{\x\)dx), (5.8) 

and let the integral kernel of the operator K he given by (2.25). Then the 2-modified 
perturbation determinant dct2(/ + -ft') admits the representation 

dct2(/ + K) = e^-Mctcd(y+(0) + r-(0)), (5.9) 

where G is defined in (4.3) and Y± are the matrix-valued Jost solutions on M± 
introduced in Definition 5.2. 
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In the remaining part of this section we will apply Theorem 5.3 to study the 
Schrodinger equation (3.1) with an integrable potential, 

V G i^(K). (5.10) 

We consider equations (1.1) and (1.2) with 

= J)' ^(-) = (nx) S)' (5-11) 

To avoid confusion we emphasize that A is an a;-independent function of the param- 
eter k in (3.1). Clearly, (1.2), with A = A{k) and R as in (5.11), is the first-order 
system corresponding to the Schrodinger equation (3.1). We note that G{A(k)) = 
{ik, —ik}, where, without loss of generality, we choose k so that Im(fc) > 0. Since 
we intend to apply Theorem 5.3, we have to make sure that Hypothesis 2.8 is sat- 
isfied. In particular, the unperturbed equation with the x-independent coefficient 
A{k) must have an exponential dichotomy on M. This is equivalent to the require- 
ment a{A{k)) n v'M = or, taking into account the choice of fc, to the inequality 
Im(A;) > 0, which we will assume to hold in the remaining part of this section. 

Let Q{k) be the spectral projection for A{k) so that cr{A{k)\^^n{Q(k))) = {ik}. 
We remark that 

x'(Q(A;)) = x(Q(fc)) = Re{ik) and ^(/d - Q{k)) = x{Id - Q{k)) = Re(-ifc). 

(5.12) 

Next, we introduce the sub-exponential weight function by 

f{x) = c, XGR+, (5.13) 

for an appropriate constant c > 1. Since the eigenvalues ±ik of A{k) are simple, 
the estimate (5.3) with $(fc,a;) — e^^^C^), a; G M, holds, and thus Hypothesis 5.1 is 
satisfied. The matrix- valued Jost solutions Y±{k,-) in the sense of Definition 5.2 
are the unique solutions of the Volterra integral equations 

/•oo 

y+(A;,a;) = e^^('=)Q(fc)- / dx' e^'^-'^'^^^''^ R{x')Y+{k,x'), x>0, 
Y_{k,x) = e^^('=)(/d - Q{k)) + r dx' e^='-''"^^^^'^R{x')Y_{k,x'), x<Q, 

J —OO 

(5.14) 

such that the matrix- valued functions Z±{k.x) = e^'^''^Y±(k,x), x G IR±, are 
bounded. Since V G i^(K), the conclusion of Theorem 5.3 holds with Y±{k,-) 
given by (5.14). 

Next, we will relate the finite-dimensional determinant in (5.9) and the classical 
Jost function (see, e.g., [6, Ch. XVII], [46, Sect. 12.1] for the latter). First, we recall 
some well-known notions from scattering theory (see [28] for a detailed bibliogra- 
phy). The Jost solutions, u±{k, •), of the Schrodinger equation (3.1) are defined as 
solutions of the Volterra integral equations 

u±{k,x) = e^'^'' - j dx' k-^sin{k{x-x'))V{x')u±{k,x'), ^^^^^ 

Im(fc) > 0, a; G M. 

The Jost function, J = J{k) is defined by 

J{k) = ^W{u_{k,-),u+{k,-)), lm(fc)>0, (5.16) 
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where W{u{x),v{x)) = u{x)v'(x) — u'{x)v{x) is the Wronskian of u and v with 
u,v € C"^(M). For fc > 0, J'(k) is the reciprocal of the transmission coefficient (see, 
e.g., formula (XVII. 1.36) in [6]). We note that J{k) is independent of x G M since 
trc2(A(A:)) =0. 

Next, we consider the factorization V = ViVr with 

Vi{x) = |y(a;)|^/^ Vr{x) = \V{x)\^/'^ exp{ia.Tg{V{x)), x gR. (5.17) 

To make the connection with (2.23), we remark that 

= o)' o)' (^-i^) 

Finally, we introduce the integral operator, L{k) on L^(R) with integral kernel 

Lik,x,x') = \Vr{x)e'''\'-'-'='Wi{x'), Im(fc) > 0, x,x' G R. (5.19) 

The following corollary of Theorem 5.3 recovers a well-known relation between 
the Fredholm determinant of the operator I + L{k) and the Jost function originally 
due to Jost and Pais [36] (sec also [44], [45], [46, Sect. 12.1.1] and, specifically, [28, 
Theorem 4.7]). More importantly, it shows that the classical Jost function coincides 
with the finite-dimensional determinant detc<i(y+(fc, 0) + Y-{k, 0)) in (5.9) obtained 
by means of the matrix- valued Jost solutions introduced in Definition 5.2. As we 
will see in Section 9, this determinant is, in fact, the Evans function associated with 
equations (1.1) and (1.2). 

Theorem 5.4. Assume that the potential of the Schrddinger equation (3.1) satisfies 

V E L^(M), and fix k G <C such that Im(fc) > 0. Let K(k) be the integral operator 
on L^(R)^ with the integral kernel given by (2.25), where $(fc,x) = e^^^''^ a; G M, 
and A{k) is defined in (5.11), and Re and Rr are defined by (5.18). Let L{k) be 
the integral operator on L^(R) with the integral kernel given by (5.19) and J^{k) be 
the Jost function of (3.1). Finally, let l±(fc, •) be the matrix-valued Jost solutions 
(5.14) on M± for (1.1) and (1.2) with A{k) and R given by (5.11). Then the 
folloiLting assertions hold: 
(i) J(fcj =detc2(r+(A;,0) + F_(A:,0)). 

{ii) The 2-modified Fredholm determinant admits the following representation: 

det2(/ + K{k)) = det2(/ + L{k)) = e®^''^J{k), (5.20) 

where 

0W = ^_^'^a;y(ar). (5.21) 

Proof. It is convenient to diagonalize A{k) and e^'^^'^^ x e K. To this effect we 
introduce the matrices 

= T(t)-'yl(fc)r(t) = (•* (5.22) 

1 / 1 {iky 



Q{k)=T{k)QT{k)-' = -{.^ [ ). (5.23) 
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We note that Q is the spectral projection for A{k) so that 

'^(^(fc)lran(Q)) = ^(^(fc) lran(Q ) = {^A:}. (5.24) 

Passing to the matrix-valued functions Y±(k,x) = T{k)~^Y±{k,x)T{k), x G R±, in 
(5.14), we obtain the integral equations 

y+(A:,x) =e^^('=)g- / dx'e(^-^')^('=V(x')5(fc)y+(fc,:r'), x>0, (5.25) 

J X 

Y_{k,x)=e''^{I-Q)+ f da;'e(^-^')^('=V(a;')S'(fc)F_(fc,a;'), a; < 0. (5.26) 

J — oo 

The property Y+{k,-) = Y+{k,-)Q and (5.25) imply that the second column of 
the 2x2 matrix Y+{k, •) is equal to zero, while Y-{k, •) = Y-{k, ■){I2 — Q) and 
(5.26) imply that the first column of Y-{k, ■) is equal to zero. Thus, in the matrix 
equations (5.25) and (5.26) we can separately consider the first and the second 
column, respectively. Let y+{k, •) denote the first column of Y+{k, •) and y-{k, ■) 
denote the second column of y_(fc,-). Passing in (5.25) and (5.26) to the C^- 
valued functions y±{k,x) = T{k)y±{k,x), x € M-t, we observe^ that they satisfy 
the following C^-valued integral equations: 

y+{k, x) = e*'="(l ikY - / dx' V{x')T{k)e^'=-'='^^^^^S{k)T{k)-^y+{k, x'), 

J X 

a; > 0, (5.27) 

y-{k,x) = e-'^'={l -ikf + f dx' V{x')T{k)e^''-=''^^<^''^ S{k)T{k)-^y-{k,x'), 

J —OO 



x<0. 



A direct calculation using (5.22) and (5.23) shows that 



T(fc)e(^-^')^('=)5(fc)T(fc)-i = 



fc-isin(fc(a;-x')) 0\ 
cos(fc(x - x')) o) ' (5.28) 

Im(fc) > 0, x,x' e R. 

Since Y±{k,-) are matrix-valued solutions of (1.2), it follows that y±{k,-) are C^- 
valued solutions of (1.2). Next, we denote by u±{k,-) the top entry of the 2x1 

vector y±{k,-). Differentiating the first components in the 2x1 vector integral 
equations in (5.27), we observe that y±{k,x) = {u±{k,x) u'j.{k,x)Y , x G IR±, 
and thus u±{k,-) are (scalar-valued) solutions of the Schrodinger equation (3.1). 
Moreover, it follows from (5.27) that u±{k, ■) satisfy (5.15), that is, u±{k, •) are the 
traditional Jost solutions of (3.1). This proves assertion (i) since 

detc2(r+(fc,o) + r_(fc,o)) = dctc2(r+(fc,o) + y_(fc,o)) 

= detc2((y+(fc,0) y-(fc,0))) = detc2(T(fc)-i(2/+(fc,0) y-(fc,0))) 

Finally, we turn to the proof of assertion (ii). A direct computation in (4.3) using 
(5.11) and the formula for Q{k) in (5.23), verifies the formula for Q{k) in (5.21). 
Using (2.25) and (5.18), and with the help of the diagonalization described in (5.22) 



■^We use T for transposition so that (a is a (2 x 1) column vector. 
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and (5.23), one computes that the matrix-valued integral kernel K{k, x, x') in (2.25) 
and the scalar- valued integral kernel L{k, x, x') in (5.19) are related by the formula 

K{k,x,x') = L{k,x,x')Q, x,x'€^, (5.30) 

where Q is the projection in (5.23). This implies det2(/ + K{k)) = det2(/ + L{k)). 
Thus, relation (5.9) and assertion (i) yield assertion [ii). □ 

While we assumed Im(fc) > throughout this section, we note that continuity of 
K{k) and L{k) with respect to fc, Im(fc) > 0, A; ^ 0, in the Hilbcrt Schmidt norm, 
permits one to extend the results of Theorem 5.4 (i), {ii) to all Im(A;) > 0, A; 7^ 0, 
by continuity. 

6. Generalized Matrix- Valued Jost Solutions 

In this section we start the discussion of a generalization of the matrix-valued 
Jost solutions introduced in Definition 3.1, and prove an extension of Theorem 3.2 
assuming that the dichotomy projection Q admits further exponential splitting. As 
we will see later on, the gcncralizcid matrix-valued Jost solutions will allow us to 
relax the exponential fall-off hypothesis on the perturbation imposed in (4.5). We 
recall from Remark 3.5 that these hypotheses are optimal, provided one merely 
assumes the existence of an exponential dichotomy Q with no further splitting. To 
simplify the exposition in this section, we will only consider the case where the 
dichotomy projection Q admits a splitting of order two and postpone the general 
case of the exponential splitting of arbitrary order until Section 7. 

Hypothesis 6.1. Assume Hypothesis 2.8. In addition, assume that the dichotomy 
projection Q can be represented in the form Q = Qi + Q2, where Qi and Q2 are 
projections which are uniformly conjugated by ^ onR such that the Bohl segment 
[h'{Qi), h{Qi)] lies strictly below [x'((32), >«'(<32)]; that is, 

x'iQi) < x(Qi) < x'(Q2) < x(Q2) < < ^(/d -Q)< >c{h - Q)- (6.1) 

We will consider the generalized dx d matrix-valued Jost solutions on K+ (cf. 
Definition 3.1), the case of M_ can be treated similarly. 

Definition 6.2. Assume R £ LI^^{RY^'^ and Hypothesis 6.1. Then d x d matrix- 
valued solutions y_|:^' and Y_^'' of the perturbed differential equation (1.2) on IR+ 
are called generalized matrix-valued Jost solutions associated with the projections 
Qi and Q2 if 

log\\Yl'\x)-^x)Q,\\c^.. ^^rn\ (a o\ 

limsup < H_^_{Qj) (6.2) 

x—*-oo ^ 

and 

Yi'\0) = Yi^\0)Qj, 3 = 1,2. (6.3) 

Wc emphasize that the additional technical requirement in (6.3) can always be 
satisfied as soon as (6.2) holds for some solutions y\^\ j — 1,2. Indeed, if the 

solutions y!^'' satisfy (6.2), then F^'^'q^, j = 1,2, are solutions of the perturbed 
differential equation satisfying both (6.2) and (6.3). 
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Remark 6.3. We note that the generahzed matrix-valued Jost solution is 
not unique. Indeed, if Y^^ is a solution satisfying (6.2) for j = 2, and 1+ is 

any solution of (1.2) whose upper Lyapunov exponent on belongs to the Bohl 
segment [>f^|_(Qi), x+((3i)], then the solution Y)^ ' +1+ again satisfies (6.2) for j = 2 
since 

log \\y!^\x) + Y+{x) - ^{x)Q'2\\c^.^ 



lim sup ■ 



X 

K2) 



^ /,. log||y|^(a:)-^(a:)Q2||c''x^ m \\ ^ 'in\ 
< max < hmsup , }t+{Qi) > < >t^{Q2)- 



(6.4) 



O 



We start with the following elementary fact. 
Lemma 6.4. Let a > and assume that < p e L^(M+). Then, 

lim e""^ / dx' e"^'p(a;') = 0. (6.5) 

Proof. Assuming mes {{x > \ p{x) > 0}) > and introducing q{x) = dx' p{x'), 
a; > 0, an integration by parts yields 



' [ dx' e"^V(a;') = q{x) - ae""^ / dx' e"^'g(a;'), a; > 0. 
Jo Jo 



(6.6) 



Since p is nonnegative and integrable on M+, lim^^oo q{x) is positive, and the 
integral dx' e"^ q{x') diverges. Using I'Hopital's rule one obtains assertion (6.5). 

□ 

The following result ensures the existence of the generalized matrix- valued Jost 
solutions on M_|_ introduced in Definition 6.2 under much weaker exponential fall-off 
assumptions on the perturbation than (3.5) and (3.6) (we recall that the require- 
ment (3.5) is optimal if we want to deal with the matrix- valued solutions introduced 
in Definition 3.1 only). 

Theorem 6.5. Assume Hypothesis 6.1 and the condition 

\\R\\cd^i&L^{R+;el^''dx) (6.7) 

with some 

/3>max{A+(g,)-xV(Q,)}. (6.8) 

J = l,2 

Then the perturbed differential equation (1.2) has generalized matrix-valued Jost 

solutions y!;^\ j = 1, 2. 

Proof. Existence and uniqueness of the generalized .Jost solution Y^^ on M_|- asso- 
ciated with the projection Qi follows from Theorem 3.2 and Remark 3.3; to check 
condition (6.3) for j = 1 we note that y!^\x) = Y^\x)Q, a; > 0, since the solution 
is unique. 

By hypothesis and the definition of the lower Bohl exponents (cf. (2.5)), for any 
£ > there exists a positive constant C{e) > such that 

\\^{x){I - Q{)^{x')-^\y.. 

< \\i^{x)Q2^{x')-^\\cixi + \\^{x){Id - Q)^{x')-^\\cix^ 
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< (eK(Q2))(x-x') ^ g(;<;(/,-Q))(c.-c.')) j^^j. _ j.'!)^ < X < x', (6.9) 
where ,fe{x) = C{e)pf^^^. Since m:'^{Q2) < >f+(^d — Q), we infer the estimate 

|j$(x)(/ - Oi)$(.t')"'||c^x. < 2e(^V(Q2))(— _ ^'D^ 0<x<x'. (6.10) 
By the definition of the upper Bohl exponent, one may also assume that 

||$(.T)gi$(a;')"'||c^x. < 2e("+('3i»(^-^'V£(|x - x'\), < x' < x. (6.11) 

Taking e > such that < 2e < (3 — {^+iQ2) — ^+iQ2)), and introducing the 
function 

p(a;) = e(^+(^2)-<(^2»^/2(,T)||i?,(.T)||c.xd, a;>0, (6.12) 

one observes that < p G i^(IR+) by assumptions (6.7) and (6.8). For brevity we 
denote fi = X+{Q2)i ol — x'+{Q2) — x+{Qi)^ and, for any r > 0, we introduce the 
integral kernel 

pru x') - e-'^C--') AM - QiMxT'Rix'), r<x<x', 

fe{x) \^X)Q^^X')-'R{X'), T<X'<X. ^ ' 

Using (6.10) and (6.11) and the inequalities x^((32) — A+(Q2) < and a > 0, one 
derives the following estimate: 

\\F;{x,x')\y.. < 2Ux)-'f^{x')\\Rix')\\c^.. 

{g(j<+(Q2)-A+(Q2))(x-a;')^ T < X < x' , 
g(^+(Qi)-A+(Q2))(x-x')^ r<x'<X, 

= 2A(a=)-ie(-;(Q^)-^+«^))^p(^') (':.(._.o " ; (6.14) 

I e ^ T < x <x, 

<2f,{x)-^p{x'), x,x'>T. (6.15) 

On the space Cb([r, oo))''^'* we define the integral operator by 

/oo 
dx'F^(x,x')Z{x'), x>T, (6.16) 

and consider the corresponding Fredholm-type integral equation 

Z{x)=Z^°\x)-iF;:Z){x), x>T, (6.17) 
where Z^°^ is defined by 

Z^°\x)=e-^^f-\x)<Pix)Q2, x>T. (6.18) 

Since /e > 1 and p G (R+ ) , the estimate in (6. 15) shows that the integral operator 
F^ is a contraction on Ch{[T,oo)Y^'^ for r > sufficiently large. Since the expo- 
nential decay of ||$ (52||c''>;'' is controlled by fi, we sec that Z^^^ G Cb([T, cxd))'*^'*. 
Thus, the integral equation (6.17) has a unique solution Z £ Cb([r, oo))'^^'' that 
can be obtained by the iteration process 

oo 

Z{x) = 2'(^)(x), Z'-^^x) = (i^;z(^-i))(x), j eN, x>T. (6.19) 

i=i 

Clearly, by uniqueness of the solution, 

Z{x) = Z{x)Q2, x>T. (6.20) 
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A straightforward computation shows that the function 

Yl^\x)=e^'-Mx)Z{x), x>T, (6.21) 

yields a solution of the perturbed differential equation in (1.2). Moreover, from 
(6.17) and (6.14) it also follows that 



(6.22) 



Mx)\\Z{x) - Z(0)(a;)||c.x. < 2eK(«^)-^+('3=))- (^|b|Ui([,,^)) 

+ 6-"==^ rfx'p(x')e"'^'^||Z||Cb([r,oo))<i>«i, X>T. 

Using p S L^(M+) and Lemma 6.4, (6.22) implies the asymptotic relation 

e—'+(Q2)-\\Yl^)^a:)-^{x)Q2\y>.. = o(l). (6.23) 

Moreover, using (6.7) and (6.8) again, choose 6 S (0,a) sufficiently small such that 
the function ps{x) = e^^p{x), a; > 0, is integrable on M+. Then (6.22) implies 

fe{x)\\Z{x) - zW(a;)||c.x. < 2e-'^e^<^Q^^-^+^Q^^^^ (\\ps\\lh[.,oo)) 

. ^ (6.24) 

+ e-("-^)^ (ia;'w(x')e("-')"' j ||^||c,([r,oo))<^x., x > t. 

Using ps G L^(M_|_) and Lemma 6.4 again, one obtains the asymptotic relation 

e-<(Q-)-\\Yl^\x)-^{x)Q2\\c^.^ = o(e-*"). (6.25) 

X^OO 

(2) 

Finally, the matrix-valued solution FJ: ^ can be uniquely extended to the interval 
[0, r) by solving the initial value problem 

yI^^'{x) = {A{x)+B{x))yI^\x), x>t, Yl^\T) = e^^UT)Z{T). (6.26) 

In addition, from (6.25) one derives (6.2), and from (6.20) it follows that Y_^\o) = 
Y^\o)Q2, completing the proof. □ 

In order to discuss uniqueness properties of the generalized d x d matrix-valued 

(2) 

Jost solution Fj: up to lower-order terms (cf. Remark 6.3), we recall that if Q is 
the (unique) exponential dichotomy projection for the unperturbed equation (1.1) 
on M, then equation (1.1) also has an exponential dichotomy on R+ with the same 
dichotomy projection Q. However, the corresponding Bohl segment on M+ may be 
strictly smaller than that on M, that is, [k'^{Q), >c+{Q)] C [V(Q), >c{Q)], see (2.10). 
By Lemma 2.13, the perturbed differential equation with ||i?||c<ixd G L^(R+) will 
also have an exponential dichotomy on IR+ with a dichotomy projection P having 
the same Bohl exponents as the unperturbed differential equation for the dichotomy 
projection Q. We emphasize that the dichotomy projection P is not unique: Only 
the subspace ran(P) is determined uniquely since it consists precisely of those initial 
data uo such that the C^-valued solutions of the initial value problem 

y'{x) = {A{x) + R{x))y{x), x>0, y{0) = yo, (6.27) 

are bounded on M_|_. Moreover, by Lemma 2.14, if Q = Qi + Q2 is an exponential 
splitting for (1.1) on R (and therefore on M+), then the dichotomy projection P for 
(1.2) on R_|_ also admits an exponential splitting P = Pi + P2 with the same Bohl 
segments, that is, [>/|_(Pj), x+(Pj)] = [>c'_^{Qj), H+{Qj)], j = 1,2. The projection 
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Pi is not uniquely determined but its range is, since ran(Pi) consists of precisely 
those; initial data yo such that the Lyapunov exponent of the C'-valued solution 
of (6.27) satisfies A+(y) < >(+{Pi) < >c'^{P2) = x'+{Q2) (cf. Remark 2.7). Thus, 
in view of Theorem 6.5 and Remark 6.3, we have proved the following uniqueness 
result: 

Corollary 6.6. Assume the hypotheses of Theorem 6.5. Then: 

(i) The generalized matrix-valued Jost solution Y_^^ associated with the projection 

Qi is unique. 

{ii If y!^^ and v!^^ are any two generalized matrix-valued Jost solutions associated 
with the projection Q2 then ran(y_[^^(0) — ^^^^(0)) C ran(Pi). 

7. Exponential Splitting of Arbitrary Order 

In this section we consider the generalized matrix- valued Jost solutions of (1.2) in 
the general case where the unperturbed equation (1.1) has an exponential splitting 

of arbitrary order on M. 

Hypothesis 7.1. Assume Hypothesis 2.8 with d > 2. Suppose, in addition, that 
for some d' , 2 < d' < d and ko, 1 < ko < d' — 1, the dichotomy projection Q 
on R for (1.1) admits an exponential splitting Q = X]j=i Qj order feg, and the 

projection {Id — Q) admits an exponential splitting I — Q = Yl'j=ko+i '^■f order 
d' — ko, where the projections Qj, j = 1, 2, . . . , d', are uniformly conjugated by i> 
on M and the corresponding disjoint Bohl segments [>c' {Qj), x{Qj)\ are ordered as 
follows: 

x'{Qj) < x{Qj) <>i'{Qj+i) < x{Qj+i), l<3<d'-l. (7.1) 

Definition 7.2. Assume R G Ll^^{R)'^^'' and Hypothesis 7.1. Then d x d matrix- 
valued solutions Y^^\ j = 1, . . . ,kQ, on M+ and y}_-'\ j = ko + 1, . . . , d' , on M_ of 
the differential equation 

Y'{x) = {A{x) + R{x))Y{x), X G M±, (7.2) 

are called generalized matrix-valued Jost solutions associated with the exponential 
splitting {Qjjf^i if 

Hmsup '°^ll^"W-^Wftllc'" < ^^(g^), (,.3) 

yI'\0) = yI'\0)Q„ j = l,2,...,fco, (7.4) 

and 

^^^yo,\\Yi^\^)-H^)QAy^^^ ^_ ^^^^^ 

x^ — 00 X 

F«(0) = yW(0)Q,-, i = fco + l,...,d'. (7.6) 
Assuming Hypothesis 7.1, and referring to Lemmas 2.13-2.14, let 

P+ = Pi + P2 + • • • + Pa-„ (7.7) 
be the exponential splitting of order ko of the dichotomy projection P+ for (1.2) on 
R+ such that 

W+{Pj), x+(P,-)] = W+{Qj), x+(Q,)], j = l,...,ko. (7.8) 
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Similarly, let 

P_ = Pfe,+i + Pk,+2 + --- + Pd' (7.9) 

be the exponential splitting of order {d' — ko) of the dichotomy projection P_ for 
(1.2) on M_ such that 

K(P,),>^-(P,-)] = K(Qj),^-(Qj)], j = fco + i,...,d'. (7.10) 

Our next result shows that the generalized matrix-valued Jost solutions exist and 
arc unique up to lower exponential order terms. 

Theorem 7.3. Assume Hypothesis 7.1 and the condition 

||i?||c<ixd G Li(R;e''l^lda;) (7.11) 

with some 

/? > max I max {X+iQj) - <(Q,)}, ^ niax {x_(Q,) - X'_{Qj)}}. (7.12) 

I l<j<ko ko+l<J<d' ) 

Then the following conclusion holds: 

{i) The perturbed equation (1.2) has generalized matrix-valued Jost solutions in the 
sense of Definition 7.2, j = 1, 2, . . . , fco, on M+ and y'^\ j = ko + 1, ■ . ■ ,d' , 
on M_ associated with the exponential splitting {Qj}'j^i. 

For any set of generalized matrix-valued Jost solutions in the sense of Definition 
7.2, the following assertions hold: 

{ii) The generalized matrix-valued Jost solutions vj^^ and Y^"^ ^ , associated with the 
projections Qi and Q41, are uniquely determined and satisfy 

ran(y|^^(0)) C ran(Pi), ran(Yl'''^(0)) C ran(Pd.)- (7-13) 

{Hi) The generalized matrix-valued Jost solutions Y^\ j = 2,...,ko, and Y^^\ 
j = ko -\- 1, . . . ,d' — 1, satisfy 

ran(yj^)(0))\{0} Cranf'^Pfc'jV ran(^PfcY j=2,...,ko, (7.14) 

ran(yi^)(0))\{0} Cranj'^Pfc'jy ran( ^ Pfc) , j = ko + 1, . . . ,d' - 1. 

(7.15) 

(w) //rj^'^ andvl^^ are any two generalized mairix-valued Jost solutions associaied 

(i) ~ (i) 

with the projections Qj for j = 2, . . . , ko and Y_ and Y_ are any two generalized 
Jost solutions associated with the projection Qj for j = fco + 1, . . . , rf' — 1, then 

ran (yf^'^ (0) - Y^'^ (0)) C ran ^ P^ V j = 2, . . . , ko, (7.16) 

^ fe=i ^ 

ran {y['^ (0) - y['^ (0)) C ran ^ P^ V j = ko + 1, . . . ,d' - 1. (7.17) 

^ k=j + l ' 

Proof. We will provide a sketch of the proof for the case of referring for details 
to the proof of Theorem 6.5. The case of IR_ is treated similarly. 
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To prove (i) and the uniqueness statement in (m), we fix j e {!,..., /co} and 
denote = A+(Qj) +£, j = 1, . . . , fco, for some £ > 0. In addition, for some r > 0, 
we introduce the integral operator F^^^ with integral kernel 



F;/(x,x') = e-''^(- 




^{x'y^Rix'), T <x<x', 

(7.18) 

'S>{x')-^R{x'), T<x'<x, 



on the Banach space Ch{[T, (x))'^^'^ . Repeating the arguments in the proof of 
Theorem 6.5 under hypothesis (7.12), one concludes that the operator F^j-' is a 
contraction on Cb([r, oo))'^^'*, provided r is large enough and e is sufficiently small. 
Thus, for each such r, the Fredholm-type integral equation 

Z«(x) = 4^)(x)-(F;/zW)(a;), x>r, (7.19) 
where Z^^ is defined by 

Z^^\x) = e-^^''^x)Qj, x>T, (7.20) 

has a unique bounded matrix-vahied solution on [t, oo) that can be obtained by 
iterations similar to (6.19). We note that if j = 1, then equation (7.19) is a Volterra 
integral equation, and there is no need to pass to a large r to ensure the contraction 
property: Indeed, the Volterra integral operator F^'^^ has zero spectral radius on 
Cb([T, oo))'^^'' for any r > and one can start the iteration process to obtain the 
solution even at r = (cf. the proof of Theorem 3.2). In particular, this shows the 
uniqueness part of (ii). If j = 2, . . . , fco, then the generalized Jost solution Y_^^ on 

M-l- is first constructed on the interval [r, oo) by Y^^\x) = e^^^^ Z^^\x), x > t, and 
then extended to [0, r) by solving the initial value problem 

{y}'^)'{x)^{A{x)+R{x))y}_'\x), a;G[0,T], yI'\t) = ef"'^ Z^^\t). (7.21) 

One verifies as in the proof of Theorem 6.5 that (7.3) and (7.4) hold. (We note the 
fact that for j = 2, . . . , /co, the generalized matrix- valued Jost solutions depend on 
r). This yields the existence of the generalized Jost solutions associated with the 
exponential splitting Q = 5^j^Li Qj '^^^ tlie dichotomy projection Q. 

Inclusion (7.13) in (ii), and assertions (Hi) and (iv) follow from (7.3), (7.4) 
and the elementary properties of Lyapunov exponents of C^-valued solutions of 
differential equations listed in Remark 2.7. Indeed, turning to the proof of (7.13) 
and (7.14), we first remark that for J = 1, . . . , ko, 

y log\\Yl'\x)\y.. ^ r \og\\Yi'\x) - <S>(x)Q,\\c^.. ^ ^^^\ 
limsup < max < limsup ,X^(Qj) > 

X — >oo X y x^OO ^ ) 

(7.22) 

does not exceed >c^{Qj) = x+(Pj) by (7.3). Thus, the Lyapunov exponent A+(y) 
of the C'-valued solution y of (1.2) with initial data y(0) e ran {y\^'' (0)) is strictly 
smaller than AY(-Pj+i) for j = 1, . . . , fco — 1 and is negative for j = ko- Applying 
(2.19) to the perturbed equation (1.2), we obtain the inclusion ran (y^-''' (0)) C 
ran ( J2k=i for j = 1, . . . , /cq- To finish the proof of (7.14) for j = 2, . . . , fco, 
let us suppose that there is a nonzero vector j/o = Y^^\o)yi, with some yi e C'', 
such that yo G ran ( J2k^i ^k) • By (7.4) we have Qjyi 7^ 0. Then the Lyapunov 
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exponent of the C^-valued solution y of the initial value problem (6.27) given by 
y{x) = Y^^\x)yi, x >0, must satisfy A+(y) < x+(Pj_i) by Remark 2.7. This leads 
to a contradiction proving (7.14). Indeed, adding and subtracting y{x) = Y_^\x)yi 
in ||$(a;)Qj2/i||c<i, we conclude from (7.3) that 



lim sup 



log{\\<S>{x)Qjyi\\cd) 



log||y[^")(x)-cl>(x)g,||c.x. 1 
<inax< limsup ,x+(Pj_i) > 

x^oo ^ ) 

< x'4Qj). (7.23) 

But this is impossible since the Lyapunov exponent of the nonzero C^-valued so- 
lution ^Qjyi of (1.1) on K_|_ with initial data Qjyi G ran(<5j) must belong to the 
Bohl segment [K'_^{Qj), K+{Qj)] (see again Remark 2.7). 

Finally, to verify assertion (iv) we add and subtract ^{x)Qj in Y^-'\x) — Yj_''\x) 
and use (7.3) to conclude that A+(?/) < X'_^_{Pj) for every C'-valued solution y of 
(1.2) with y(0) e ran (0) - y}_^^ (0)) . This proves inclusion (7.16). □ 

Remetrk 7.4. In addition to assertions {ii)-{iv) in Theorem 7.3, the generalized 

matrix- valued Jost solutions , constructed in the existence part (i) of Theorem 
7.3 by means of the Fredholm-typc integral equations (7.19) and their analogs 
for have the following property: If n e N, is the sequence of truncated 
perturbations as in (4.14), and Y^\ denote the generalized matrix- valued .lost 
solutions associated with the perturbation i?„ (see the proof of Theorem 4.1), then 

lim ri^'^(O) = yp\0). (7.24) 

To see this one follows the course of the proof of Theorem 6.5. First one establishes 

continuity of the map R Z^^^ from L-'^(R+; e^'^dx) to Cb([T, oo))''^'*, defined in a 
neighborhood U{R) of R in the space L^(IR+; e^^dx). Subsequently, one shows that 
the convergence of the sequence {iin}neN from this neighborhood to i? as n — > oo 
in the L^{^+; e'^^c/a;)-topology yields (7.24). O 

Next, given any set of generalized matrix-valued Jost solutions, wc introduce 
the following Evans determinant (the terminology is related to the Evans function, 
which is further discussed in Section 9). 

Definition 7.5. For a given set of generalized matrix- valued Jost solutions Y^\ 

(i) 

j = l,...,ko, on IR+ and Y_ ', j = ko + 1 ■ ■ ■ ,d , on M_, the Evans determinant, 
D, is defined by 

D = detc4y+ +y-), where y+ = ^r[^)(0), 3^- = ^ ^^'^(O). (7.25) 

The following important and purely algebraic result shows that although the 

generalized matrix-valued Jost solutions are not unique, the Evans determinant D 
is uniquely determined by equations (1.1) and (1.2). 
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Lemma 7.6. Assume Hypothesis 7.1. Then the determinant D in (7.25) is in- 
dependent of the choice of the generalized matrix-valued Jost solutions Y^^ , j = 
l,...,ko, and Y^^\ j = ko + 1 . . . ,d' . 

Proof. Using the notation (7.7) and (7.9), let P+ and P_ bo the dichotomy projec- 
tions and {-Pjl^li and {Pj}'^^ko+i ^^'^ exponential splittings for (1.2) on and 
M_, having the same Bohl segments as the splitting {Qj}j^i for (1.1). Without 
loss of generality, we will assume that the subspaces ran(Pj), j = 1, . . . ,d' , satisfy 
ran(Pi)-i- ran(P2)+ • • • + raii(Pc(/) = (otherwise the determinant D is equal to 
zero for any choice of the system of the generalized Jost solutions). We note that 
ran(3^±) C ran(P±) by assertions (ii) and (Hi) in Theorem 7.3, and y+ = y+Q 
and y_ = y-{Id - Q) by (7.4) and (7.6). Thus, treating the matrices y± in (7.25) 
as operators 3^+ : ran(<5) — > ran(P+) and : ran(Jd — Q) ^ ran(P_), we ob- 
serve that the block-operator y = y+ + is diagonal in the following direct sum 
decomposition: 

y = (^^ : = ran(Q)+ ran(/rf - Q) ^ = ran(P+)+ ran(P_). (7.26) 

Assertions {ii) and {Hi) of Theorem 7.3 then also yield that the matrix of the 
operator 

y+ : ran((5) = ran(Qi)-i- • • • -i-ran((5fcj 

. . (7.27) 
^ ran(P+) = ran(Pi)-F ■■■+ ran(Pfe„ ) 

is upper- triangular, while the matrix of the operator 

y- : ran(/d - Q) = ran{Qko+i)+ i-ran(Qd/) 

(7.28) 

ran(P+) = ran(Pfco+i)-| |-ran(Pd/) 

is lower-triangular. Moreover, from assertion {iv) in Theorem 7.3 it also follows 
that the diagonal blocks of the operators y± with respect to the decompositions 
(7.27) and (7.28) are independent of the choice of the generalized Jost solutions, 
completing the proof. □ 

Remark 7.7. In fact, the proof of Lemma 7.6 only uses properties (7.4) and (7.6), 
and properties (7.13) (7.15) of the matrix-valued solutions of (1.2). O 

At this point we are ready to prove the principal result of this paper. 

Theorem 7.8. Assume Hypothesis 7.1 and the condition 

WRWcdxd e L\R;e'^^''^dx) (7.29) 

for some 

p > max I max {X+{Qj) - >c'+{Qj)}, ^ max {x_(Q,) - A'_(Q,)}}. (7.30) 

Let Yj_-'\ j = 1,2, . . . , ko, and Y^-'\ j = ko + 1,2, . . . , d', be any system of matrix- 
valued generalized Jost solutions. Let K he the integral operator on L^(M)'' whose 
integral kernel is given by (2.25). Then the 2-modified perturbation determinant 
det2(/ -|- Jl') admits the representation 

det2{I-hK) = e®detc4y++y-), (7-31) 
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where Q is defined in (4.3), y± are defined in (7.25), and D = detc£i(3^+ + 3^-) is 
the Evans determinant {of. (7.25)). 

Proof. Let Y^^l^, j = 1, . . . , fco, and vi"']?, j = ko + 1, . . . ,d' , be the generalized 

matrix- valued Jost solutions on R_|_, respectively, on R_, constructed in the ex- 
istence part (i), {a) of the proof of Theorem 7.3 for a given r > sufficiently 
large. Here, the subscript F is added to remind the reader that these solutions 
are obtained from solutions of the Fredholm-type integral equations (7.19) and 
their analogs for M^. In general we note that the solutions Y^^\p depend on r, 

and they are not the solutions given a priori in the formulation of Theorem 
7.8. Introduce the truncated perturbations n e N, by formula (4.14), and 
denote by Y^^]^ p, j = 1, . . . ,ko, and Y^-']^ p, j = fco + 1, • ■ • , d', the correspond- 
ing generalized matrix- valued Jost solutions of the truncated perturbed equation 
y'{x) = {A{x) + Rn{x))y{x), obtained by solving the Fredholm-type integral equa- 
tions (7.19) with R replaced by Rn (see the proof of Theorem 7.3 (i)). Since R^ is 
compactly supported, there exist unique matrix-valued Jost solutions in the sense 
of Definition 3.1 of the truncated perturbed equation on M_|_ and on R_; these will 
be denoted by and y_,„. Using the solutions F±,„ and the projections Qj, 
define matrix-valued solutions of the truncated perturbed equation by 

y+lvi^^) = Y+,n{^)Qj, J = l,...,ko, X>0, 
y-,l,vi^) = y-,n{x)Q3, j = ko + l,...,d', X<0. 

Here, the subscript V is added to remind the reader that the solutions y±,„ are 
obtained from the Volterra integral equations (4.26). Since the support of i?„ is 
compact, equations (4.26) imply that the solutions Y^l^ y constitute a system of 
generalized matrix- valued Jost solutions in the sense of Definition 7.2. In addition 
to (7.25), introduce the notations 

/CO ^' 

y+,F =E>^IS(0)' ^-.^ = E ^-!f(0), (7.33) 



i=i 



j=fco-l-l 



y+,n,F =E41f(0), y-,n,F = ^ ^-'IfW, (7.34) 
■'-^ j=ko + l 

ko 

■'-^ j=ko + l 

Using (4.27) and the definition of Y^]^ y in (7.32) we obtain 

fcp 

y+^,,y = J2 y+,n{0)Qj = i + ,n(0)Q = Y+,n{0), (7.36) 

i=i 

and similarly, 

y_,n,v = Y-,n{0). (7.37) 
Using Remark 4.3 this results in 

det2(/ + K) = e® lim detc^ {y+,n,v + y-,n,v) ■ (7.38) 
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(i) (7) 

Since both Y^'^ y and Y^'^ p constitute systems of generalized matrix- valued Jost 
solutions, Lemma 7.6 yields 

detcd {y+,n,v 

By the continuity argument discussed in Remark 7.4 (cf. (7.24)), we get 

lim detci {y+,n,F + y-,n,F) = detc<i {y+,F + y-,F) , (7.40) 
and next, again by Lemma 7.6, 

detc. {y+,F + y-,F) = detcd {y+ +y-). (7.41) 

Combining (7.38)-(7.41) yields (7.31), completing the proof. □ 

Next, wc will give an extension of Lemma 7.6 and Theorem 7.8 for the case 
where the matrix-valued solutions ^Qj, j = l,...,d', of (1.1) in Definition 7.2 
are replaced by arbitrary solutions of (1.1) with certain asymptotic properties. We 
assume Hypothesis 7.1 and fix the splitting {Pj}j^i for (1.2) described in (7.7) 
and (7.9). Let Nj, j = 1, . . . ,d', be any matrices in C'^^'' satisfying the following 
conditions: 

(i) Nj=NjQj=QjNj, j = l,...,d'. 

(7.42) 

{ii) detc<i(iV) ^ 0, where N = 22^i- 

3 = 1 

Let $j for j = 1, . . . , fco, respectively, j = ko + 1. . . . , d' , denote the matrix- valued 
solutions of the unperturbed equation (1.1) on K+. respectively, M_. satisfying 
the initial conditions ^j(0) = Nj, j = l,...,d'. By («) in (7.42), the Lyapunov 
exponents of the C^-valued solutions of (1.1) given by the columns of belong 
to the Bohl segments corresponding to the projections Qj, j = l,...,d'. Next, 

assume that Y^^K j = 1, 2, . . . , fco, on R+, and y!:^\ j = ko + 1,2, . . . ,d' , on 
are any given d x d matrix- valued solutions of the differential equation Y'{x) = 
{A{x) + R{x))Y{x), X G M±, satisfying the following properties (cf. Definition 7.2): 

\og\\Yi^\x) - ^Ax)\\cd>^i , , , , , 

lim sup ^" + ^ ^ ^-^-^ — < >i'+{Qj), (7.43) 

a:— »oo X 

Y!f\Q) = Y!f\Q)Qj, j = 1,2,..., fco, (7.44) 

and 

x^-oo x 

yi^) (0) = y['^ (0)Q„ j - fco + 1, ■ ■ • , d'. (7.46) 
Define the following matrix N (cf. Definition 7.5): 

fco d' 

7V = 5;++3;_, where 5^+ = 2?+^(0)> 5^- = ^-^W- (7-47) 

We remark that if Nj = Qj in (7.42), then detc£i(A^) = 1 and Yj^^ are generalized 
matrix- valued Jost solutions i'^ the sense of Definition 7.2, and N = y+ + y-, 
in the notations of (7.25). 
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Lemma 7.9. Assume Hypothesis 7.1. In addition suppose that Nj satisfy (7.42), 
that $j satisfy $j(0) = Nj, and that Y^^^ satisfy conditions (7.43) -(7.46) for j = 

Finally, assume that N is defined as in (7.47). Then the ratio of the 
determinants dct^d (N) and det^d (N) is independent of the choice of the solutions 
with the initial data satisfying (7.42) and the choice of the solutions j = 

1, . . . ,d'. Thus, one obtains for the Evans determinant D in Definition 7.5, 

detc<i(iV) 

Proof Using (ii) in (7.42), we denote yI^\x) = yI^\x)N-'^, x € R±, for j = 
I,. ..,d . Wc claim that Y^ arc generalized matrix- valued Jest solutions in the 
sense of Definition 7.2. Indeed, since Qj = NjN~^ due to Nj = QjN, the equalities 

yI^\x) - ^x)Qj = yI^\x)N-^ - ^x)NjN-^ = {yI^\x) - ^j{x))N-^ (7.49) 

show that condition (7.3) follows from (7.43). Also, since N~^ and Qj commute, 

condition (7.44) implies (7.4), proving the claim. By Lemma 7.6, the Evans de- 
terminant D constructed using the solutions Y^^^ as described in Definition 7.5, is 
independent of the choice of Nj and 1± , and hence equals 



(feo d' s 

J2yI^\o)n-^+ y['\o)n-^] 



(7.50) 



= detcd(7VAr-i), 

proving (7.48). □ 

Corollary 7.10. Under the assumptions imposed in Theorem 7.8 and Lemma 7.9, 
the 2-modified perturbation determinant det2(/ + K) admits the representation 

where 6 is defined in (4.3), and N and N are defined in (7.42) and (7.47). 

Remark 7.11. Theorem 7.3 can be viewed as a further development of the cele- 
brated Levinson theorem (see, e.g., [13, Chap.l] and the bibliography cited therein). 
The Levinson theorem for the asymptotically diagonal perturbed equation (1.2) on 
R+ deals with the situation when A{x) = dia,g{Xj{x)}j^i, x G R+, and ||ii||c<ix<i G 
L^(]R+) and asserts the existence of C*- valued solutions of (1.2) satisfying the as- 
ymptotic relation 



Jo 



yj{x) = {ej+o{l))exp Xj{s)ds, j = l,...,d, (7.52) 



where ej, j = 1, . . . ,d, are the vectors of the standard basis in C^. The underlying 
assumption in the traditional Levinson theorem (i.e., for (7.52) to be valid) is that 
the following alternative holds: For each pair of integers j,k e {1, . . . , d}, j ^ fc, 
either 

liminf 7j,/j(0,a;) = — oo and sup Ij^k{x' ,x) < oo, (7.53) 

x^oo x>x'>0 
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or 



a;>a;'>0 



inf k(x, x') > —00, 



(7.54) 



where we denote Ij^k{x',x) = J^, ds Re(Aj(s) — Xk{s)), < x' < x. To compare 
our results and the traditional Levinson theorem, we first note that, unlike the 
assumptions in the Levinson theorem, we do not assume that A is diagonal. In the 
class of diagonal unperturbed equations (1.1) our assumptions on the unperturbed 
diagonal system are more special (since we are interested only in bounded solutions 
on IR+, and assume the exponential dichotomy on M), but in turn we derive more 
conclusions in addition to the asymptotic relation (7.52), see assertions {ii)-{iv) 
in Theorem 7.3. If, in addition, we assume that the diagonal unperturbed system 
has d disjoint Bohl segments, then the above mentioned alternative holds; thus 
our assumptions are stronger than that in the Levinson theorem in this particular 
case. However, generally, our hypotheses are more flexible since we can group the 
diagonal elements of A related to the same Bohl segment, thus avoiding conditions 
on all pairs j,k as in the traditional Levinson theorem. Our assumptions (7.11) 
on the perturbation are stronger than those in the Levinson theorem, but in turn 
(cf. (7.52)), we conclude that the columns of the generalized matrix- valued Jost 
solutions approximate the reference solutions up to terms o{e~^^) as x — > cxd for a 
positive 5 (cf. (7.3) and (6.25)). Finally, the generalized matrix- valued Jost solutions 
yield formula (7.31) for the perturbation determinant, a feature that is not discussed 
in the context of the traditional Levinson theorem. O 

8. Autonomous Perturbed Equations 

In this section we treat the case of an autonomous unperturbed equation (1.1), 
that is, we consider the differential equations 



on M, where A G C'^^'' is a constant matrix and R G ^^^^(M)''^'' (cf. Examples 2.3, 
2.4), and 2.6). Our objective is to show that the conclusions of Theorems 7.3 and 
7.8 hold for (8.1) and (8.2) under much weaker assumptions on the perturbation 
R than the exponential decay imposed in (7.11) and (7.12). Specifically, we will 
merely assume that ||i?||c£ix£i G L^(M) when the eigenvalues of A are semi-simple, 
or ||i?||cdx£i decays polynomially when A has non-diagonal Jordan blocks. Under 
these weaker assumptions, the asymptotic behavior of the generalized matrix-valued 
Jost solutions will no longer be measured in terms of exponential weight factors as 
indicated in Definition 7.2. Instead, we will define the generalized matrix-valued 
Jost solutions of (8.2) by means of solutions of a certain mixed system of Volterra- 
and Fredholm-type integral equations in the spirit of Definition 5.2. 

Hypothesis 8.1. Assume thatd >2,A€ C^^"^, and a{A)r\iM. = 0. LetQj denote 
the spectral projection of A corresponding to the spectral subset of the eigenvalues 
of A having equal real parts denoted by Xj, j = l,...,d', where 2 < d' < d. In 
addition, assume that for some feg, 1 < ko < d' — 1, the inequalities ki < ■ ■ ■ < >ffej, < 
< Xko+i < ■ ■ ■ < ^d' hold, so that Q = Qj is the dichotomy projection for 

(8.1) on M and on ]R±. 



y'{x) = Ay{x), a; G R, 

y\x) = {A + R{x))y{x), a; G M, 



(8.1) 
(8.2) 
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Passing to an appropriate coordinate system, we may assume that the matrix A 
is in Jordan normal form, and thus for each j = 1, . . . ,d' the operator ^|ran(Q,) is 
represented by a direct sum of diagonal matrices vl and/or matrices 1/1 + J, where 
v e (^{A), Re(i^) = Kj, and J is the matrix of an appropriate size with I's above the 
main diagonal and all other entries equal to zero. We introduce nij £ N such that 
nij + 1 is equal to the maximal size of the non-diagonal Jordan blocks of ^|ran(Qj)j 
noting that rrij > 1 if such blocks exists, and setting nij = if all Jordan blocks are 
diagonal (i.e., if all eigenvalues u of A with Re(z/) = xj are semi-simple). Clearly, 
the matrix exponent e'^''"'"'''^' can be computed explicitly using the Jordan blocks 
of vl|j.a,n(Qj) (see, e.g., [14, Example 1.2.5]), so that with this notation one has the 
estimate 

||e"^^l'-«^)|| < ce^^^(l + |a;|)™^ j = l,...,d', xgR, (8.3) 

for some constant c > 0. Finally, we introduce m = maxi<j<d/{mj}. 

Next, we fix T > 0, and consider the following "mixed" system of Volterra- 
and Fredholm-type integral equations (cf. (5.6) with the weight function f{x) = 
(1 + la;!)™^, xGR, and (7.19) with /x^ = xj): 

Z^/\x) - (1 + |x|)-"^e^(^-"^'Qj 

= - / da;'e(^-^')(^-"^)(l + |a;|)-'"^ ( (1 + WD""' R{x')Z^^\x') 

+ £ dx' e(--')(^— ^)(l + |ar|)-™' + Wl)""' R{x')Z^^\x'), 

j = l,...,ko, x>T, (8.4) 

Z^J\x) - (1 + |a;|)-™^e^(^-''^-)Qj 

= r dx'e(^-^')(^-'^-')(l + |x|)-'"^Y^QfcVl + |a;'|)'"^i?(a;')^i'^(a;') 

- f "da;'e(^-^')(^-"^Hl + kl)""^( Qk){l + Wir^ R{x')Z^_!\x'), 

j = fco + 1, . . . , rf', x< -T. (8.5) 

Here, we set X^^Z^ Qk — ^ when j = 1 in (8.4) and Qk — ^ when j — d! in 

(8.5); thus, the first and the last equations in (8.4) and (8.5) are of Volterra-type, 
and the remaining equations are of Fredholm-type. 

Definition 8.2. Assume Hypothesis 8.1 and R £ L\^J^'^'^'^ . Then d x d matrix- 
valued solutions y'^^ , j = 1, . . . , fco, on R_|_ and y'^^ , j = /cq + 1, . . . , d', on R_ of 
the differential equation 

Y'{x) = {A^R{x))Y{x), a;GR±, (8.6) 

are called generalized matrix-valued Jost solutions associated with the exponential 
splitting {Qj}jLi if 

Yi^\x) = {l + \x\r^e''^''Z^^\x), i = l,...,fco, x>T, (8.7) 
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= (1 + |a;|)™^e''^="zi^'^(x), j ^ ko + 1, . . . , d' , x < -t, (8.8) 

where Z^_^\ j = 1, . . . , A:o, are bounded solutions of equation (8.4) on [r, oo) and 
j = fco + 1, . . . are bounded solutions of equation (8.5) on {—oo,t]. The 
solutions Y^^ and Y^^^ are extended to [0, r) and to (— r, 0], respectively, by solving 
the initial value problem for equation (8.2) with initial data Y^\±t) = (1 + 
T)"^i e^'^i'^ Z^\±t) for the corresponding values of j = 1, . . . , d'. 

Wc emphasize that the; generalized matrix-valued Jost solutions in the sense of 
Definition 8.2 are not unique (unless d' = 2, see Definition 5.2) and depend on the 
choice of r in (8.4) and (8.5). We will continue to use Definition 7.5 of the Evans 
determinant, but the generalized matrix- valued Jost solutions in this definition will 
be understood in the sense of Definition 8.2. We recall the exponential splitting 
{Pj}j^i in (7.7) and (7.9), for which the Bohl and Lyapunov exponents for the 
perturbed equation (8.2) are also equal to hj, j = 1, . . . ,d' . 

Theorem 8.3. Assume Hypothesis 8.1 and the condition 

\\R\\c-ixi e i^(K; (1 + |a;|)2™ dx), (8.9) 

where m+\, m,> 1, is the maximal size of the non-diagonal Jordan blocks of A 
[if such blocks exist), or m = if all Jordan blocks of A are diagonal. Then the 
following conclusions hold: 

Assertions {i)-{iii) of Theorem, 7.3 hold, where the generalized matrix-valued Jost 
solutions are understood in the sense of Definition 8.2. 
In addition: 

(iv) The generalized matrix-valued Jost solutions satisfy 

e—o-\\Yl^\x)-e-^Qj\y.. = o(l), j = l,...,ko, 

... ^ ™ (8.10) 

X— > — oo 

{v) The 2-modified perturbation determinant det2{I -\- K) admits the representation 
det2{I + K) = e^detc4y+ + y-), {8.11) 

where 

POO 1*0 

e= / dxtrc4QR{x)]- dxtic4{Id-Q)R{x)], (8.12) 

Jo J -oo 

y± are defined as in (7.25), and D = detcd(3^+ + y-) is the Evans determinant 
{of (7.25)). 

Proof. We will sketch the proof for the case of IR+ referring for details to the proofs 
of Theorem 7.3, Lemma 7.6, and Theorem 7.8. The proof for the case of R_ is 
similar. We fix j € {1, . . . , fco}. Starting the proof of assertion («), let FJ^'^^ denote 
the integral operator on Cb([r, oo))'^^''' defined by the right-hand side of (8.4) (cf. 
(7.18)). 

We claim that under assumption (8.9) the following assertions hold: 
(a) FJ^'^^ is a contraction on Cb([T, oo))'^^'^ for t sufficiently large. 
(6) There exists a function g e Cb([T, oo)) so that limj^^oo g{x) = and 

\\{FlfZ){x)\\c^., < g{x){l + |a;|)-™^ ||Z||c,([.,oo))<*x. (8.13) 
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for all X e [r, oo) and Z £ Cb([T, oo)f''. 

To prove the claims (a) and (6) we use (8.3) and bound the norm of the first integral 
in (8.4) by the expression 

C|l4''||c.([r,cx>))<^x4l + kl)-'"^- 

X / dx' ^e("'=-"^)("-"')(l + |x-x'|)""=(l + |x'|)"^||E(a;')||c<ix. 
"^^ fc=j 

/•OO 

<c|lz(f'|lc,([.,oo))-x. / dx'(l + |a;'|)2'"^||i?(a;')||c''x., (8.14) 

since, due to x < x', for fc > j + 1 in the sum above, we have Xk — > and 
exp((>!r/j — >fj)(x — x'))(l + |x — x'l)™* < c, while for k = j we have (1 + |x — x'|)"'j < 
c(l + |x' I )'"■' . Similarly, the norm of the second integral in (8.4) is dominated by 
the expression 

c|l4'^||c.([r,oo))-<^(l + k|)-"^ 
fX J-1 

X / rfx' ^e("'=-"^-'(^-^')(l + |x-x'|)'"'=(l + |x'|)"^||i?(x')||c^xd 
"'^ fe=i 

< c||Z^^^||crj,([^,oo))<ix<i 
/■X J-1 

X dx' ^e("'=-"^-'(^-^')/2(l + |x'|)2™^-||i?(x')||c''x<i, (8.15) 
"'^ fc=i 

since cxp((>i:fc — >fj)(x — x')/2)[l + |x — x'!]™*" < c due to Xk — >ij < and x — x' > 0. 
Now condition (8.9) and Lemma 6.4 yield claims (a) and (6). 

Using claim (a), assertion (i) follows. Also, the uniqueness part in assertion (ii) 
holds since (8.4) is a Voltcrra integral equation for j = 1. Assertions (7.13), and 
{Hi) in Theorem 7.3, and assertion (iv) in Theorem 8.3 follow from claim (6) similar 
to the proof of Theorem 7.3. For instance, to prove (Hi), we note first that claim 
{b) implies the estimate 



\yI^\x) -e^'^QiWcdxd <cg{x)e''i'', x > r, where lim ^(x) = 0. (8.16) 

X — >QO 



If y is a nonzero C'-valued solution of (8.2) with y{0) e ran(F["'^ (0)) then, using 
(8.16), its Lyapunov exponent satisfies the inequality 

> / ^ f, logiirpVx) -e^'^g,-||c<ixd 1 w X 

A+(y) < max<^ limsup ^'11 + ^^ 'im^ — K < ^.^^ = A+(P,+i), 

j = 2,...,fco-l, (8.17) 
A+(t/)<0, j = ko. (8.18) 

Applying Remark 2.7 to (8.2) we conclude that y{0) G ran ( ^^^-^ Pfc) for j = 
1, . . . , fco. However, for j = 2, . . . , ko, the assumption y{0) = S ran ( X]i=i ^k) 
leads to a contradiction similar to the proof of Theorem 7.3 {Hi). Indeed, using the 
same notation as in that proof, we have 

We^^QmWc < c\\e^^Qj-Yl^\x)\\c<ix^ + \\y{x)\\c<i < C5(x)e-^-^ + ce-^-^^ (8.19) 
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where we used (8.16), the assumption y{0) e ran ( Ffe) , and Remark 2.7 
applied to (8.2). It follows that e~''^^\\e''^Qjyi\\c<i — > as a; — > oo, in contradiction 
with the explicit formula for g^"^'""*'')' in terms of the Jordan blocks (cf. [14, 
Example 1.2.5]). This proves assertion {in). Assertion {iv) follows from (8.16). 

It remains to prove assertion (v). For this purpose we consider the truncated 
perturbed equation, y'{x) = {A + Rn{x))y{x), x € M., with i?„ as in (4.14). For 
a sufficiently large t > 0, we define the operator F^'^!'.^ on Cb([r, cxd))'^^'^ by the 
right-hand side of equation (8.4) with R replaced by Rn- Since i?„ — > i? in 

L^(R; (l+bl)^'" dxY""^, estimates similar to (8.14) and (8.15) show that F;>.J„ ^ 

^ ' ra— ►oo 

FJ^'j' in operator norm on the Banach space Cb([T, oo))''^''. It follows that 

YiiM^) - r£k0), (8.20) 

(i) 

where Y^!^ p are the generalized matrix-valued Jost solutions corresponding to the 

truncated perturbed equation, and Y^^ = Y^p are the generalized matrix-valued 
Jost solutions of (8.2) constructed in part (i) of the current proof. The rest of the 
proof is similar to the arguments in the proof of Theorem 7.8. Indeed, using the 
notations introduced in that proof, to establish formula (8.11), we need to show 
that det2{I + K) = e®detc<i(y+,F + 3^-,f)- The latter equality follows from (7.38), 
(7.39), and (7.40). As before, (7.38) follows from Remark 4.3, and wc already know 
from (8.20) that (7.40) holds. Thus, it remains to prove (7.39). For this we will 
take advantage of the fact that the support of Rn is compact (and thus we may 
drop the improper integral in (8.4) for sufficiently large a; > 0). Indeed, we may 
assume that n > r. It follows from (4.26) and (7.32) that Yl^\ y{x) = e^^Qj for 
x> n. On the other hand, for a; > n we also have (cf. (8.4)) 

y+lA^) = e^^Qj + J^x' e^^-^'^^(j2Q'^)Ri^')YiiA^^ (8-21) 
Thus, using (8.3) for x > n and a sufficiently small £ > 0, we obtain the estimate 

dx'^e-(---')(l + |x-x'|)"- 
fe=i 

X \\R{x')\\c^.4Yl^lpix')\\c^.. 

<c/ dx' Ve^'"'"^'^'^"'^'^ max e'^'^+^^^^-^'Vl la; - a;'|)"'= 
~ Jr fr[ x>x'eml ^ ' 

X \\R{x')\\c^.4yI^Ip{x')\\c^.. 
< ce'^'^i-^> with c = c(n, r), (8.22) 

since Xk < Xj for k = l,...,j — 1. Using Remark 2.7 we therefore obtain the 
inclusion ran (y|^i^^(0) - Fj^i,y(0)) C ran ( ECl Pfc), J = 2, . . . , k^. Now (7.39) 

follows as in the proof of Lemma 7.6 since Y_^^!^^p and Yj^^!^^ possess all properties 
{ii)-{iv) in Theorem 7.3 (see Remark 7.7). This concludes the proof of assertion 
{v) in Theorem 8.3. □ 
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Remark 8.4. We emphasize that ahhough the generahzed matrix-valued Jost 
solutions in the sense of Definition 8.2 are not uniquely determined and depend on 
the choice of r in (8.4) and (8.5), formula (8.11) shows that the Evans determinant 
D is independent of the choice of the solutions and hence independent of r. O 

Remark 8.5. We note that Theorem 8.3 {iv) also follows from the Levinson theo- 
rem (cf. Remark 7.11) for asymptotically constant coefficients under the additional 
assumption that the spcctrTim of A is simple. In this context we refer to [13, The- 
orem 1.8.1] and also to [13, Sect. 1.10] for additional results on the asymptotic 
behavior of solutions of (1.2) with A of Jordan block-type. O 

9. The Evans Function 

In this section we relate the Evans determinant D introduced in Definition 7.5 
to the Evans function^ E. First, we recall one of many equivalent definitions of 
the Evans function available in the literature (see, e.g., [1, 5, 35]), namely, the 
definition using exponential dichotomies, see [57, Definition 4.1]. Consider a family 
of diff'erential equations 

y'{x)=B{z,x)y{x), a; e K, (9.1) 

parameterized by a complex spectral parameter z 6 O, where J7 C C is open and 
simply connected. It is assumed that the locally integrable C''^''-valued functions 
B{z, ■) depend on 2 e O analytically. Since in this paper we are not concerned^ 
with function theoretic issues related to the Evans function such as its analytic 
continuations, etc., we will fix a value zq G O in what follows, and will give the 
definition of E = E{zo) for this Zq G ft. Accordingly, we will suppress the z- 
dependence in the notations for the differential equations above, and consider just 
one differential equation y'{x) = B{x)y{x). 

Hypothesis 9.1. Assume that the differential equation y'{x) = B{x)y{x), a; G R, 
has an exponentially hounded propagator on M and, in addition, an exponential 

dichotomy P+ on IR+ and an exponential dichotomy P_ on R_ . Moreover, assume 
that the ranks of the projections P+ and P- are equal, and denote the common 
value by n so that 

n = dim(ran(P+)) = dim(ran(P_)) and d — n = dim(ker(P+)) = dim(ker(P_)). 

(9.2) 

Definition 9.2. Assume Hypothesis 9.1. Choose any vector basis {ei, . . . ,e„} of 
the subspace ran(P+), and any vector basis {e„+i, . . . , e^} of the subspace ker(P_). 
The Evans function, E, is then defined as the determinant of the dxd matrix with 
columns ek, k = 1, . . . ,d. 

Clearly, E = ii and only if the subspaccs ran(P+) and ker(P_) have a nonzero 

intersection. C'-valued solutions on M+ with initial data y+'(0) = ek, k = 
1, . . . ,n, form a basis in the subspace of all bounded solutions of the differential 

equation y'{x) = B{x)y{x), x G R+ . Similarly, C'-valued solutions y^^^ on R_ 

with initial data y^^\Q) = e^, k = n + 1, . . . ,d, form a basis in the subspace of 
all solutions bounded on 1R_. Thus, E = if and only if the differential equation 
has a C'-valued solution bounded on M; see the Introduction. We emphasize, that 



^See, however, Remark 9.5 below. 



EVANS FUNCTIONS, JOST FUNCTIONS, AND FREDHOLM DETERMINANTS 



45 



the Evans function E just defined is not unique, and depends on the choice of the 
vectors e^, fc = 1, . . . , d. Moreover, this definition does not assume any perturbation 
structure in the differential equation, and, in particular, can be used to define the 
Evans functions for both (1.1) and (1.2). 

Returning to the principal theme of the current paper, wc will set B(x) = A{x) + 
R{x), X gR, with the matrix-valued functions A and R as in (1.1) and (1.2). Next, 
we intend to show that, in fact, D = E, where D is the Evans determinant in 
Definition 7.5, and E is the Evans function defined using a special choice of the 
vectors e^, fc = 1, . . . , d, given by the initial data of the generalized matrix- valued 
Jost solutions of (1.2). 

For this, we assume in accordance with Hypothesis 2.8, that the unperturbed 
dxd matrix differential equation y'{x) = A{x)y{x), x gR, with A G Ll^^(R)'^^'', has 
an exponentially bounded propagator and an exponential dichotomy Q on R. As- 
sume, in addition (see Hypothesis 7.1), that the unperturbed equation (1.1) has an 
exponential splitting of order d' , 2 < d' < d, so that = ran((5i)+ • • • + ran((5d')! 
and that the Bohl segments corresponding to the projections Qj are disjoint. We 

also have Q = X^^li Qj and I — Q = X^^^^g+i Qj for some 1 < fco < d'. Sup- 
pose that a. d X d matrix-valued function R is such that ||i?||c£ixd G _L'^(R). Un- 
der these assumptions, due to Lemma 2.13 and its analog for we conclude 
that the perturbed equation y'{x) = {A{x) + R{x))y{x), x G IR±, has an expo- 
nential dichotomy P-f on R+ and an exponential dichotomy P_ on R_ such that 
dim(ran(Q)) = dim(ran(P+)) and dim(ker(0)) = dim(ker(P_)). Thus, (9.2) holds 
if we let 

fcp 

n = dim(ran((5)) = ^ dim(ran((5j)). (9.3) 

The same conclusions hold provided the unperturbed equation (1.1) is autonomous, 
and we assume Hypothesis 8.1. We summarize our assumptions as follows. 

Hypothesis 9.3. For the unperturbed equation (1.1), respectively, (8.1) assume 
Hypotheses 2.8 and 7.1, respectively, Hypothesis 8.1. Also, assume that ||i?||c<ixd G 
L-^(IR). In addition, assume that the perturbed equation (1.2), respectively, (8.2) 
has a system of generalized matrix-valued Jost solutions y\^\ j = 1, . . . , ko, on R_^_ 

and Y^^\ j = ko + 1, . ■ ■ ,d' , on IR_ in the sense of Definition 7.2, respectively, in 
the sense of Definition 8.2. 

Our final result then reads as follows. 

Theorem 9.4. Assum,e Hypotheses 9.3. Define the m,atrices y± and the Evans 
determinant D as indicated in Definition 7.5. In addition, with n as in (9.3), let Cfe, 
fc = 1, . . . , n, denote the nonzero columns of the matrix y+ and e^, fc = n-|- 1, . . . , d, 
denote the nonzero columns of the matrix . Letting B{x) = A{x) -\- R{x), x G K, 
and using the vectors ek, k = 1, . . . ,d, define the Evans function E for the perturbed 
equation (1.2) as indicated in Definition 9.2. Then 

E = D. (9.4) 

Moreover, if. in addition, the perturbation R satisfies the assumptions in Theo- 
rem 7.8, respectively. Theorem 8.3, then the 2-modified perturbation determinant 
det2(/-|-^^) admits the representation 

det2{I + K) = e®E, (9.5) 
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where G is defined in (4.3), respectively, in (8.12). 

Proof. We claim that n nonzero columns yk{x) of the d x d matrix y^{x) = 
Ylj=i X > 0, provide a basis in the space of all bounded C*- valued so- 

lutions of the perturbed equation (1.2) on K_|_, while d — n nonzero columns yk{x) 
of the matrix y-{x) = Yl'j=ka+i^-^ ^ — yield a basis in the space of all 
bounded C'-valued solutions of the perturbed equation y'{x) = {A{x) + R{x))y{x) 
on M_. Indeed, the claim follows from (7.3) and (7.5) (or, for the autonomous 
unperturbed equation (1.1), from assertion (if) in Theorem 8.3) and the fact that 
the nonzero columns of the matrix ^{x)Q, respectively, <J?(a;)(/d — Q), provide a 
basis in the space of all bounded solutions of the unperturbed equation (1.1) on 
M+, respectively, K_. 

As a result, letting e^ = 2/^(0), k = 1,. . . ,d, we selected a basis ei, . . . , e„ in 
the subspace ran(P+), and a basis e„+i, . . . ,erf in the subspace ker(P_) needed for 
the definition of the Evans function E as above. In particular, this proves that 
the Evans determinant D equals the Evans function E defined with this choice of 
bases. The result in Theorem 7.8 (or, for the autonomous unperturbed equation 
(1.1), Theorem 8.3(f)) now states that, under additional assumptions on R (see 
(7.11) and (7.12) or (8.9)), the 2-modified perturbation determinant det2(/ + K) 
admits the representation (9.5). □ 

Remeirk 9.5. Assume that the matrix-valued functions A{z, •) and R{z, •) in (1.2) 
analytically depend on a spectral parameter ^ G f2 C C. Moreover, assume that 
the equation y'{x) = A{x, z)y{x), x G M., has an exponential dichotomy on M 
for each z G fl. Then the dichotomy projection Q = Q{z) for (1.1) on M. is an 
analytic C''^''-valucd function on O as well. This follows, for instance, from the 
discussion in Remark 2.11 (alternatively, see [53] or [57, p. 995]). Thus, Q{z) in 
(4.3) is an analytic function for z € fl. Moreover, the corresponding function 
K{z) is also analytic as a function with values in the space of Hilbcrt Schmidt 
operators equipped with the ,B(Z/^(R)'^)-norm. Hence, the 2-modified perturbation 
determinant det2(/ -|- K{z)) is also analytic for z G Cl. As a final result, it follows 
from formula (7.31) in Theorem 7.8 (or from formula (8.11) in Theorem 8.3) that 
the Evans determinant D{z) also depends on z G analytically. In view of the 
discussion in this section we just proved that the Evans function E{z) is an analytic 
function of the spectral parameter z G ft (cf. also [1] for an entirely different proof 
of this result) . O 

Remark 9.6. The assertions in Lemma 7.9 and Corollary 7.10 (which are more 
general than Lemma 7.6 and Theorem 7.8) can be used in the context of this section 
as follows. Let ep\ i = 1,. . . , dim(ran((3j)), denote an arbitrary basis in ran(Qj), 
j = 1, . . . , d', and let N denote the dxd matrix with the columns ep^ . For brevity, 
we enumerate the columns as e^. A; = 1, . . . , d. We note that {ek}k=i is a basis in C*. 
Passing from the standard basis in C'' (i.e., (1, 0, . . . , 0), . . . , (0, . . . , 0, 1)) to this new 
basis, we observe that for each j = 1, . . . ,d' , the matrix ((Z^fc^)^^^ ^.^^ corresponding 
to the operator Qj in the basis {cfe}J?^j^ is a diagonal matrix with all entries being 
equal to zero except the entries qfj^ = 1 for ^ = rij^i + 1, . . . , nj. Here we denote 
no = and rij = ^^^^^ dim(ran((5fe)) for j = l,...,d'. Next, we introduce the 
C'-valued solutions of the unperturbed equation (1.1) by yk{x) = ^{x)ek, a; e M, 
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k = 1, . . . ,d. Let yfc(x), x e K, k = 1, . . . ,d, be (not necessarily unique) solutions 
of the perturbed equation (1.2) that satisiy the asymptotic conditions 

Vkix) = yk{x) + o{e<^'^^^), k = l,...,n, 

x—>oo (Q (i\ 

Mx) = =yfc(a;) + o(e--(^<^-«)^), k = n + l,...,d, 

X — ^ — oo 

where n = dim(ran((5)). Now consider two Evans functions, E^y^y and E^y^y 
for the perturbed and unperturbed equations, (1.2) and (1.1), respectively. Here, 

^{vk} = dctcd(A^), and = detcd(iV) with the matrix N defined in (7.47) and 

constructed using the matrices Nj = NQj = QjN, j = 1, . . . ,d' . The representa- 
tion 

det2(/ + i^) = e®|^^ (9.7) 
now follows from Corollary 7.10. O 

Appendix A. Operators with Semi-Separable Kernels 

In this appendix we recall some of the results derived in [28]. Let di,d2 € N, 
d = di + d2, and let fj and gj, j = 1,2, be given matrix- valued functions on M 
satisfying 

/j e L2(M)'^^''^ gjeL^iRp""^, i = l,2. (A.l) 
On the space L'^{R)'^ we consider the Hilbert-Schmidt integral operator defined by 

{Ku){x) = J dx' K{x,x')u{x'), MGi^(M)'', (A.2) 

R 

where K(-,-) is a semi-separable dx d matrix- valued integral kernel defined by 

Kix,x') = {^;\'\^f^ ^f;, x,x',R. (A.3) 

[f2{x)g2{x'), x<x', 

In addition, we introduce the d x d matrix-valued integral kernel 

H{x, x') = /i(.T),gi(.T') - /2(x)52(x'), (A.4) 
and the corresponding Volterra integral equations 

/■CSO 

h{x) = h{x)-\ dx' H{x,x')A{x'), (A.5) 

Mx) = f2{x) + r dx' H{x, x')h{x'). (A.6) 



We note that equations (A.5) and (A.6) have a unique pair of solutions satisfying 
fj e L'^{RY^'^\ j = 1,2. Finally, we introduce the d x d matrix- valued function B 

by 

Bix)=('^'^f{\^f, ''^^HffX XGR. (A.7) 

V"-92G^)./lG^) -92{x)f2{x)J ' 

Theorem A.l ([28]). Assume (A.l). Then: 

(i) The first-order dx d matrix differential equation 

U'{x) = B{x)U{x), x€M, (A.8) 
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permits an explicit particular solution given by the formula 

^^^^"1^ !^ dx' g,{x')h{x') I,,-r.^dx'g,{x')h{x'))' x G M. (A.yj 

(m) The modified Fredholm determinant det2 (/ — K) has the following represen- 
tation: 



(A.10) 



det2(/ — -ft') = detcn(f/(— oo)) exp dxti:c'^{fi{x)gi{x)) 

= detc"(t/(oo))exp dxtic^{f2{x)g2{x))j . 
By Liouville's formula (cf., e.g., [33, Theorem IV.1.2]) one infers for any x,Xo € 

M, 

detcci(C/(x)) = dctcd(C/(.xo))cxp dx' tYc4B{x'))^ , (A.U) 
and thus, for any xq S R, one has from Theorem A.l {ii), 

det2(/ - K) = detc<i(C/(oo)) exp ^ J dxtrcd{f2{x)g2{x))j 

= detcd(C/(.xo))exp ( / dxtTc''igiix)fi{x) - g2{x)f2{x)) 



X exp dxty:c4h{x)g2{x))^ . {A.U) 
This results in the representation 

det2(/ - K) = detc<i(t^(a;o))exp ( / dxtrcd{f2{x)g2{x)) 



+ f dxtTc4fi{x)gi{x))\ (A.13) 

Jxo / 

= det^d {U (xq)) exp ^tr^d dx K{x — 0,x) 

+ J dxK{x + 0,x)^y (A.14) 

We note that if if is a trace class operator with a continuous integral kernel, then 
trc<i(-B(a;)) = 0, and applying [12, p. 1086-87] the Fredholm determinant can be 
represented in the simpler form 

det(/ -K)= detcd{U{xo)) = detcd{U{-oo)) = detcd(J7(oo)). (A.15) 

Appendix B. The Proof of Lemmas 2.13 and 2.14 

Proof of Lemma 2.13. Using rcscaling (sec Remark 2.2), wc will assume without 
loss of generality that x+{Q) = —>c'_^_{Id — Q)- We choose k so that < k < 
—H+{Q). We claim that: 

(i) dim A'' = dim(ran(Q)), 
and 

(ii) for a constant C(x) > solutions y of (1.2) satisfy: 

||y(a;)||c. < C(x)e— ||y(a;')||c<^, x>x'>0, if y(0) G N, (B.l) 
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||y(:c)||c. <C(;<)e''(^-^')||y(a;')||c-, > x > 0, if 2/(0) G ker(g). (B.2) 

As soon as these claims are proved, the assertions in the lemma follow. Indeed, 
since C'' = A'' © ker((5) by (i) and (B.1)-(B.2), the exponential boundedness of 
the propagator of (1.1) implies that P is an exponential dichotomy for (1.2) by 
the proof of Lemma IV. 3. 2 in [10]. (We note that in the proof of this lemma only 
the exponential boundedness of the propagator has been used.) Since — x can be 
chosen arbitrary close to M:-^.{Q), it follows from (B.l) and the definition of x+{P) 
(cf. (2.5)), that x+{P) < >cj^{Q). Considering (1.1) as an -perturbation of (1.2), 
we obtain the opposite inequality; the second equality in (2.31) is proved similarly 
using (B.2). 

Turning to the proof of the claim, let r > 0, and denote by F'^ the integral 
operator on Cb([T, oo))'' with the integral kernel F{x,x')R{x') for x,x' > r, where 
F{x, x') is defined in (2.30). By a direct estimate, inequalities (2.12) and condition 
lli^llcixd e L^(R+) imply that F'^ is a contraction on Cb{[T,oo))'^ for r sufficiently 
large. Thus, the formula 

y{x) = ^x)q + {F-'yXx), q G ran(Q), x e [r, oo), (B.3) 

gives a one-to-one correspondence between the solutions y, y{0) G A^, of (1.2) 
bounded on K+, and the solutions $g, q € kei{Q), of (1.1) bounded on In 
particular, (i) in the claim above holds. Also, due to (B.3), the bounded solutions 
y of (1.2) admit the representation (cf. [10, p. 181]) 

y{x) = ^{x)Q^x')-^y{x') + {F=''y){x'), x>x'>t. (B.4) 

Taking e > such that a = — >?+((5) — e > and applying the exponential 
dichotomy estimates in the right-hand side of (B.4), we arrive at the integral in- 
equality 

/.oo 

\\y{x)\\c<i<a{e)e-'^^^-^'^\y{x')\\c^+c{e) / dse-"l--^l||i?(s)||c.x.||y(s)||c.. 

J x' 

(B.5) 

Thus, the function u{x) = e^^||y(x)||cd, x > x' > satisfies the inequality 

roc 

u{x) < a(£)e-"(^-^')+^^||y(a;')||c<* + c(e) / dse-"l^-"l+^(^-")||i?(s)||c<ixdw(s). 

Jx' 

(B.6) 

Claim (B.l) can now be proved using arguments in [10, Sect. III.2]. Indeed, on the 
space Cb([T, oo)) we consider the integral operator T with integral kernel 

T(a;,a;') = c(£)e-«l^-^'l+''(^-^')||ii(a;')||c<i>«i, x>x'>t. (B.7) 

Using equation (2.15) in [10, Sect. III. 2] one then shows that ||T|| is bounded by 
the expression 

px+l 

c(£)[l-F(l-e-("-^))-^-F(l-e-("+^))-i] sup / ds\\R{s)\\cdx4. (B.8) 

a:G [t.oo) J x 

Thus, the operator T is a contraction for r = t{x,s) > suficiently large. Passing 
to equality in the inequality (B.6), solving the resulting integral equation for u 
by applying (/ — T)~^, and using [10, Lemma III. 2.1], we derive (B.l) from the 
estimate 

«(a;)<a(£)(l-||T||)-W||2/(a;')||c<^, x > x' > r. (B.9) 
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Turning to the proof of (B.2), we note that any solution y of (1.2) on M+ with 
y(0) e kei{Q) solves the equation 



yi'x) = ^x){Id - QMxT'yix') + I ds ^x)Q^s)-' R{s)yis) 

(B.IO) 







Jx 



ds $(a;)(7d - Q)^s)~'^R{s)y{s), 0<x<x', 
which yields the integral inequality 

\\yix)\\c^<a{e)e'''----'^\\y{x')\\c^+c{e) T dse"(--^)||i?(s)||c.x.||y(s)||c<i 



+ c{s) [ rfse-«(--«)||ii(s)||c.x.||y(s)||c<. 

Jo 

0<x<x'. (B.ll) 

Arguments similar to the proof of (B.l) conclude the proof of Lemma 2.13. □ 

Proof of Lemma 2.14- We will use the rescaling from Remark 2.2 as follows: For 
each j = 1, . . . , d' — 1 we fix /ij G (>f+((5j), >f^((5j+i)), and for j — d' we fix /i^' > 
^+{Qd)- Then, for each j = l,...,d', the rescaled unperturbed equation (2.11) 
with fj, = fij has an exponential dichotomy on R+ with dichotomy projection Qj = 
X]fe=i Qk- By Lemma 2.13, the rescaled perturbed equation, y'{x) = {A{x)+R{x) — 
I2jld)y{x), X € M., has an exponential dichotomy Vj, where Pj is the projection 
parallel to ker(Qj) on the subspace Nj consisting of the values y{0) of the bounded 
solutions of the rescaled perturbed equation on R_|_ . We note that Vd' = Id and 
set Vo = 0. Letting Pj = Vj — Vj-i, j = I,. ■ . ,d', finishes the proof of Lemma 
2.14. □ 
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